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1.  PROBLEM  DEFINITION 


1 . 1  SCOPE 

1.1.1  Background 

This  report  documents  the  second  phase  in  a  study  of  minimizing  environ¬ 
mental  effects  on  radar  sensors.  This  phase  involves  the  development  of  en¬ 
vironment-adaptive  radar  techniques  that  will  materially  enhance  the  performance 
of  a  ground-based  tactical  radar. 

The  baseline  radar  is  intended  to  handle  a  realistic  airborne  threat  that 
is  much  more  demanding  than  the  threat  that  is  handled  by  current  radars  like 
the  new  G.E.  AN/TPS-59  radars.  The  major  features  of  this  threat  are: 

1.  Large  search  volume — 360°  azimuth  by  20  km  altitude  by  200  km  range. 

2.  Multiple  target  tracking — >  50  at  2  Hz  update  rate. 

3.  Small,  high-speed  targets — <  1  sq  m  RCS,  1000  mps  speed. 

4.  Heavy  clutter —  -15  dB  terrain,  4  mm/hr  rain,  100,000  lb  chaff  cloud. 

A  conventional  design  will  not  meet  the  needs  of  this  realistic  environment. 
However,  as  long  as  the  design  is  flexible  enough  to  allow  adaptive  adjust¬ 
ments  on  the  basis  of  eventual  measurements  of  the  environment,  degraded  per¬ 
formance  in  the  most  extreme  environments  may  be  acceptable.  We  have  previously 
developed  baseline  designs  for  meeting  this  environment  without  the  aid  of  real¬ 
time  adaptivity.  In  this  report,  however,  we  develop  the  tools  for  a  fully- 
adaptive  radar  system  in  which  the  signal  processor  and  waveform  generator  adapt 
in  real-time  to  sensed  changes  in  the  environment. 

The  design  procedures  that  are  used  here  go  far  beyond  the  procedures  com¬ 
monly  found  in  the  literature.  The  paramount  problem  in  our  demanding  environ¬ 
ment  is  that  unacceptable  blind-speed  intervals  exist  when  the  usual  design 
approaches  are  followed.  Fundamentally,  the  product  of  the  range  interval  and 
range-rate  interval  to  be  handled  is  so  large  that  no  single  pulse  burst  is 

adequate.  In  general  such  problems  exist  whenever  the  cited  product  exceeds 

7  2  8  2 

cA/4  (which  is  10  m  /s  at  2.25  GHz).  Since  the  product  is  2  x  10  m  /s  for 

our  severe  environment  (R  =  200  km,  v  =  1000  mps)  the  blind-speed  prob- 
lem  overpowers  nearly  every  other  design  consideration.  Consequently  the  com¬ 
mon  design  procedures  that  are  found  in  the  literature  are  totally  inadequate 
because  they  do  not  provide  a  solution  to  the  blind-speed  problem. 


In  our  baseline  designs,  however,  we  found  that  the  problems  could  be 
solved  without  reverting  to  UHF.  By  considering  how  waveform  agility  can 
benefit  the  radar  system  we  have  developed  a  design  approach  that  greatly  re¬ 
lieves  the  remaining  problems.  This  approach  involves  a  recognition  that  the 
environment  at  short  ranges  includes  ground  clutter  and  a  need  for  rapid  track 
updates  while  the  environment  at  long  ranges  does  not.  The  consequence  of  the 
approach  is  that  a  schedule  of  waveforms  that  handles  all  ranges  can  be  de¬ 
veloped  that  is  much  more  efficient  than  the  usual  approach  which  involves  a 
single  waveform. 

1.1.2  Outline  of  Report 

In  this  report  we  first  review  the  two  baseline  systems  which  bound  the 
reasonable  alternatives  for  the  defined  environment.  The  first  system  is  based 
on  a  set  of  mechanically  rotating  arrays  with  electronic  scan  in  elevation 
only  (1-D  phased  array).  The  second  system  is  based  on  a  set  of  stationary 
arrays  with  electronic  scan  in  both  azimuth  and  elevation  (2-D  phased  array). 

The  first  system  is  very  limited  in  its  tracking  capabilities  because  of  a 
lack  of  flexibility,  but  it  is  retained  for  consideration  because  it  is  so 
much  cheaper  than  the  second  system.  In  the  process  of  reviewing  our  past  de¬ 
sign  we  also  develop  the  tools  for  an  automated  design  procedure.  The  resulting 
computer  programs  are  detailed  (and  listed  in  the  Appendix).  These  programs  are 
useful  not  only  for  easily  generating  new  designs  for  new  conditions  but  also 
for  inclusion  In  a  real-time  adaptive  system. 

Before  proceeding  to  the  adaptive  radar,  however,  we  consider  two  new  de¬ 
sign  concepts  that  appear  to  have  potential  for  further  performance  improvements. 
One  involves  the  simultaneous  use  of  multiple  carrier  frequencies  while  the  other 
involves  waveforms  where  blind-range  zones  as  well  as  blind-speed  zones  must  be 
resolved.  The  problems  with  each  technique  are  revealed  in  the  report  of  results. 

Next,  the  design  for  a  real-time  adaptive  radar  is  developed.  Automatic  al¬ 
gorithms  for  determining  the  appropriate  PRF  schedule  are  developed,  using  the 
results  of  the  previously-developed  design  tools  for  determining  the  unique 
design  parameters.  Performance  improvements  to  be  gained  from  the  use  of  such 
adaptivity  are  then  illustrated  with  a  simple  example. 

Finally,  the  plan  for  the  remaining  two  tasks  of  the  program  are  reviewed. 

It  is  concluded  that  the  next  task,  which  is  a  simulation  task,  will  result  in 
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a  final  radar  design  and  a  specification  of  the  capabilities  of  an  associated 
environmental  sensor.  The  final  task  will  then  involve  the  design  of  this  sen¬ 
sor  and  its  integration  with  the  radar  system. 

1.2  PERFORMANCE  SPECIFICATIONS 

The  baseline  systems  are  designed  to  meet  the  requirements  for  a  tacti¬ 
cal,  ground-based,  mobile,  air  surveillance  radar  in  a  mountainous  European 
environment.  These  requirements  are  reviewed  here,  for  use  in  examples  that 
appear  throughout  the  report. 

The  parameters  used  for  the  baseline  designs  are  listed  in  Table  1-1.  The 
physical  constraints  are  set  by  the  need  for  mobility.  Up  to  four  antenna 
faces  of  the  specified  size  are  allowed.  The  two-way  antenna  efficiency  is 
consistent  with  -35  dB  sidelobes.  The  surveillance  parameters  are  determined 
by  the  envisioned  tactical  environment  which  could  include  attacks  from  all 
sides.  The  cited  search  and  track  update  rates  are  specified  for  the  short 
and  medium  ranges.  Lower  rates  are  permissible  at  far  ranges.  The  target  pa¬ 
rameters  are  based  on  small  tactical  aircraft,  but  a  missile  threat  should  be 
evaluated  as  well.  The  clutter  environment  is  based  on  severe  but  not  un¬ 
likely  conditions.  The  cited  terrain  clutter  is  typical  of  dense  forests. 

The  rain  parameters  are  typical  of  heavy  storms  but  not  tornados.  The  chaff 
is  an  amount  that  could  be  carried  by  a  large  transport  aircraft. 

When  a  baseline  design  approach  involving  two  search  zones  is  considered, 
there  is  a  need  for  separately  specifying  the  parameters  in  each  zone.  In 
the  near-range  zone  we  have  clutter  consisting  of  the  combination  of  ground 
clutter,  rain,  and  chaff.  The  total  clutter  therefore  has  a  combined  Doppler 
spread  of  20  mps.  The  search  frame  time  is,  as  in  Table  1-1,  2  sec.  In  the 
far-range  zone  we  have  clutter  consisting  only  of  rain  and  chaff  clutter, 
since  the  ground  clutter  is  beyond  the  horizon  at  far  ranges.  The  total 
Doppler  spread  of  the  clutter  is  therefore  only  5  mps.  The  far-range  zone 
extends  from  the  radar  horizon  (typically  80  km)  to  the  200  km  maximum  range 
in  Table  1-1.  The  search  frame  time  at  these  long  ranges  can  be  as  long  as 
8  to  10  seconds,  rather  than  2  seconds.  These  sundry  parameters  form  the  basis 
for  the  designs  that  follow.  It  is  this  natural  separation  of  the  environment 
Into  two  range  zones  that  allow  us  to  develop  designs  at  higher  carrier  frequen¬ 
cies  than  has  been  the  common  practice  in  the  past. 
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Table  1-1.  Baseline  Specifications 


Equipm ent  Constraints : 

antenna  size 
transmitted  power 
antenna  aperture  efficiency 
system  dissipative  losses 
noise  figure 

Survei  1  lance  Parameters: 

azimuth  coverage 
range  coverage 

altitude  coverage 
max  elevation  angle 

search  frame  time 
track  update  rate 
t  racking  accuracy 

max  no.  of  simultaneously 
tracked  targets 
single-look  detection 
probability 

no.  of  false  alarms  per  frame 

Target  Parameters 

target  KCS 
fluctuation  model 
max  target  velocity 
max  target  acceleration 

Clutter  Parameters^ 

backscatter  coef  of  terrain 

max  rainfall  rate 

max  horizontal  rain  speed 

chaff  cloud  weight 

chaff  cloud  size 

chaff  cloud  speed 


2  m  high  by  3  m  wide 

100  kw  prime,  33  kw  avg.  radiated 

0.55 

10  dB 

3  dB 


360° 

5  km  to  200  km  in  search, 

1  km  to  200  km  in  track 
50  m  to  20  km 

30°  in  search, 

75°  in  track 

2  sec 
2  Hz 

10  m  in  range, 

0.2°  in  az/el  angle 

500 

0.90 

5 


1  sq  m 
Rayleigh 
1000  mps 
5  g 

-15  dB 

4  ram/ hr 

20  mps 

100,000  lbs 

10  km  x  10  km  x  1  Ian 

20  mps  max 
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2.  BASELINE  DESIGNS  FOR  AN  ADVERSE  ENVIRONMENT 

* 

The  baseline  designs  developed  in  our  previous  report  were  unconven¬ 
tional  in  the  sense  that  they  made  efficient  use  of  the  prior  knowledge  that 
the  range  extent  of  the  ground  clutter  is  limited.  The  approach  led  to  much 
better  performance  than  could  be  achieved  by  a  design  that  failed  to  recognize 
the  environmental  phenomena  and  used  a  single  waveform  for  all  ranges.  The 
details  of  the  baseline  design  approach  are  reviewed  again  here,  then  recently 
developed  procedures  for  automating  the  waveform  design  are  described,  and  ex¬ 
emplary  baseline  designs  are  repeated  to  illustrate  the  new  design  procedures. 
This  work  provides  the  groundwork  for  our  subsequent  development  of  real-time 
waveform  adaptivity  for  an  environment-adapLive  radar. 

2.1  BASELINE  DESIGN  APPROACH 

2.1.1  The  Blind-Zone  Problem  and  Its  Solution 

As  identified  in  the  earlier  stages  of  this  study,  the  paramount  problem 
for  a  radar  in  the  defined  environment  is  to  avoid  interference  from  ground, 
rain,  and  chaff  clutter  while  searching  over  a  broad  range  of  ranges  and 
range  rates.  It  was  concluded  that  the  uniform  pulse  burst  is  the  best  choice 
for  a  waveform  because  of  its  superior  clutter  suppression  capabilities.  How¬ 
ever,  numerous  blind  zones  exist  when  just  a  single  uniform  pulse  burst  is 
used.  The  number  of  such  zones  is 

4  R  v 

„  _  _ max  max 

blind  cA 

where  the  radar  must  handle  ranges  from  0  to  R  ,  range-rates  from  0  to  v  , 

e  max  &  max 

and  operates  at  a  wavelength  A.  Thus  the  number  of  blind  zones  exceeds  20  for 
our  assumed  environment  (R  =  200  km,  v  =  1000  mps)  when  the  carrier  fre- 
quency  is  above  2.25  GHz.  A  lowering  of  the  PRF  shifts  the  character  of  the 
zones  into  blind  speeds,  while  a  raising  of  the  PRF  shifts  the  character  of 
the  zones  into  blind  ranges.  The  blind  zones  are  illustrated  in  Figure  2-1 
for  the  case  where  the  20  zones  are  factored  into  10  blind  speed  zones  by  2 
blind  range  zones.  The  waveform  design  problem  is  then  to  find  a  sequence 

*J.  E.  Howard,  et  al3  Baseline  Designs  for  Mobile  Tactinal  Badary  MARK 
Resources  Report  179-4,  January  1980. 
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of  different  PRFs  that  have  non-overlapping  blind  zones,  but  with  sufficient 
efficiency  that  the  allowed  frame  times  and  available  power  are  not  exceeded 
because  of  the  extra  sets  of  PRFs. 

A  straightforward  approach  to  solving  the  blind-zone  problem  would  be  to 
first  choose  the  highest  PRF  that  avoids  blind  range  intervals,  and  to  next 
find  the  set  of  lower  PRFs  that,  when  taken  together,  make  all  regions  of  in¬ 
terest  free  of  blind  speed  intervals  as  well.  The  approach  follows  the  basic 

* 

procedures  outlined  by  Rihaczek.  However,  the  entire  procedure  can  be  sim¬ 
plified  into  an  optimization  based  on  just  two  parameters  summarizing  the  en- 
vironment  and  radar.  The  optimization  takes  into  account  the  limited  dwell 
times  for  a  given  beamwidth  and  update  rate  as  well  as  the  integer  properties 
of  the  pulse  bursts.  The  optimization  objective  is  to  minimize  the  allowable 
wavelength  for  a  given  search  time,  or  to  minimize  the  search  time  for  a  given 
wavelength. 

We  previously  found  that  we  could  systematically  develop  designs  for  the 
defined  problem  by  applying  this  two-parameter  optimization  directly.  How¬ 
ever,  we  also  found  that  the  optimization  resulted  in  designs  with  relatively 
low  carrier  frequencies  (below  L-band) ,  even  though  the  optimization  is  de¬ 
signed  to  minimize  wavelength.  Since  the  potential  benefits  of  a  higher  fre¬ 
quency  are  so  strong,  we  critically  studied  the  source  of  the  problem  and  found 
that  the  culprit  is  the  use  of  a  single  set  of  bursts  for  the  entire  range/ 
range-rate  space  to  be  searched.  When  we  later  separated  the  search  space  into 
two  distinct  range  zones  we  found  that  we  could  develop  acceptable  designs  at 
S-band  frequencies.  In  essense,  the  separate  design  of  waveforms  that  are  cus¬ 
tomized  to  each  zone  results  in  a  combined  system  that  requires  less  dwell 
time  per  search  direction,  for  a  1 ixed  set  of  conditions.  Higher  carrier  fre¬ 
quencies  then  become  admissible  even  though  the  corresponding  beamwidths  are 
narrower  so  that  the  available  dwell  time  per  search  direction  i6  shorter. 


A.  W.  Rihaczek,  "A  Systematic  Approach  to  Blind-Speed  Elimination,"  IEEE 
Trans. -AES,  Nov.  1973,  pp  940-947. 

G.  W.  Lank,  Radar  Design  for  Blind-Speed  Elimination  for  Search  and 
Track,  MARK  Resources  Report  179-3,  January  1980. 
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In  review,  the  baseline  approach  is  centered  on  a  division  of  the  search 
mode  into  two  range  zones  in  which  the  environment  is  distinctly  different. 

The  near-range  zone  is  defined  as  the  set  of  ranges  for  which  the  ground  is 
not  yet  beyond  the  horizon.  As  illustrated  in  Figure  2-2,  the  size  of  this 
zone  depends  on  the  height  of  the  radar  above  the  ground.  A  typical  compro¬ 
mise  range  extent  for  the  ground  clutter  for  a  mobile  tactical  radar  is  80  km, 
which  is  used  in  several  examples  in  this  report.  This  zero-Doppler  ground 
clutter  when  combined  with  rain  or  chaff  clutter  produces  a  relatively  wide 
total  clutter  spread  as  illustrated  in  Figure  2-3  and  2-4,  depending  on  the 
wind  conditions.  Since  the  ratio  of  PRF  to  Doppler  spread  must  exceed  two  in 
order  to  solve  the  blind-speed  problem  even  under  ideal  conditions,  a  wide 
spread  implies  the  need  for  high  PRFs.  A  high  PRF,  in  turn,  implies  a  short 
ambiguous  range.  However,  as  long  as  the  ambiguous  range  exceeds  the  maximum 
range  extent  of  the  ground  clutter  such  a  high  PRF  is  acceptable  for  detec¬ 
ting  targets  in  the  near-range  zone.  The  far-ixmae  zone  is  the  set  of  ranges 
from  the  end  of  the  near-range  zone  to  the  maximum  detection  range  of  interest. 
This  zone  is  free  of  ground  clutter  so  that  the  total  clutter  spread  is  much 
narrower  than  that  for  the  near-range  zone,  which  is  fortunate  because  only 
low  PRFs  are  allowed  at  the  long  ranges  involved.  Moreover,  the  search  up¬ 
date  time  (frame  time)  for  long  ranges  need  not  be  nearly  as  short  as  that  for 
near  ranges.  As  a  result  the  blind-speed  resolution  can  be  achieved  in  this 
zone  with  relative  ease. 

At  first  it  may  appear  that  the  two-zone  approach  consumes  time  rather 
than  saving  time  because  the  waveforms  for  each  zone  are  transmitted  in  series. 
In  actual  fact,  however,  the  increase  in  overall  efficiency  caused  by  custo¬ 
mized  PRFs  is  so  high  that  time  is  actually  saved.  For  example,  the  updating 
of  the  long-range  zone  can  be  slowed  down  so  much  that  only  a  small  fraction 
of  the  radar's  time  resources  is  required  for  the  long-range  waveforms.  And 
because  the  clutter  spread  is  so  narrow,'  relatively  few  PRFs  a.re  required  to 
avoid  Mind  speeds.  By  contrast,  the  near-range  zone  requires  rapid  updates 
but  the  raLio  of  clutter  spread  to  PRF  can  be  kept  low  because  the  PRFs  can 
be  much  higher  than  is  allowed  in  the  long-range  zone.  Therefore  the  number 
of  bursts  required  for  blind -speed  resolution  is  again  low.  The  final  result 
is  that,  overall,  the  time  available  for  tracking  is  maximized  while  fully 
•solving  the  search  problem,  and  the  allowable  wavelength  is  minimized. 
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Figure  2-3.  Typical  Doppler  Spread  of  Clutter  in  a 
Heavy  Storm 


envelope 


I 


i 


Doppler  Frequency  (i  mps) 


Figure  2-4.  Assumed  Model  for  Doppler  Spread 
of  Clutter 


The  baseline  design  approach  therefore  involves  the  reopt imi /at  Ion  of 
the  radar  waveforms  for  -each  new  design,  with  a  separate  set  of  optimized  wave¬ 
forms  for  each  of  the  two  range  zones.  Furthermore,  the  approach  could  be 
extended  by  similarly  dividing  the  azimuth  search  region  Into  separate  azimuth 
zones,  and  again  providing  different  waveforms  for  each  different  zone  t;o  as 
to  minimize  the  total  dwell  required  in  the  search  mode.  While  the  latter  ex¬ 
tension  was  not  a  part  of  our  baseline  design  for  a  non-adaptive  radar,  it  is 
a  part  of  our  approaches  to  be  described  later  for  an  adaptive  radar.  The 
critical  steps  in  all  of  our  radar  designs  for  the  mobile  tactical  radar  arc 
therefore  the  repeated  applications  of  uavefom  optimization  procedures  for 
blind-speed  resolution. 

2.1.2  Waveform  Optimization  for  Blind- Speed  Resolu t i on 

Since  the  waveform  optimization  must  be  used  repeatedly  in  considering 
various  design  alternatives,  especially  when  considering  an  envi ronment  adapt ive 
radar,  it  is  important  that  the  optimization  be  reduced  to  a  systematic  design 
tool.  Our  previously-reported  initial  attempts  at  a  simple  procedure  tesulted 
in  a  pair  of  design  graphs  that  could  be  used  to  roughly  estimate  the  required 
number  of  bursts  and  the  allowable  carrier  frequency  for  ihe  optimum  solution 
to  the  blind-speed  problem  in  a  specified  zone.  Since  that  time  we  have  refined 
the  graphical  procedures  to  make  them  more  useful  in  practice  and  we  have  de¬ 
veloped  new  computer  programs  that  allow  one  to  bypass  the  graphical  procedures 
and  take  advantage  of  computer-aided  design  procedures.  The  graphical  proce¬ 
dures  are  described  here  in  order  to  introduce  the  objectives  and  accomplish¬ 
ments  of  the  waveform  optimization. 

The  waveform  optimization  can  be  performed  graphically  as  follows.  First, 
two  fundamental  input  parameters  must  be  defined  from  the  environment: 
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where 

k^  “  antenna  beanwidth  factor  (which  is  0.8S  for  the  baseline  antenna) 

k  “  filter  broadening  factor  (which  is  1  to  1.7  depending  on  Doppler 
"  sidelobe  levels  required) 

v  -  maximum  target  speed  of  interest  (which  is  1000  mps  for  the 
maX  baseline 

w  “  antenna  width  (3m  baseline) 

a 

Tg  ■  total  available  dwell  time  in  one  search  frame  (less  than 
2  seconds  baseline  for  the  near-range  zone) 

0  “  angular  extent  to  be  searched  in  time  T  (2tt  baseline) 

.  8  8 
AR  ,  -  actual  range  rate  extent  of  clutter  to  be  suppressed  (20  mps 
baseline) 

The  optimization  that  minimizes  wavelength  (or,  equivalently,  minimizes  dwell 
time)  results  in  two  parameters  which  uniquely  define  the  solution: 

K,  the  ratio  of  maximum  target  Doppler  to  the  maximum  PRF  in  the  set  of  PRFs 
p,  the  number  oi  separate  bursts  required. 

This  optimum  can  be  found  directly  from  Figures2-5  and  2-6.  The  solution  can 
then  be  translated  into  radar  parameters  as  follows.  The  basic  (maximum)  PRF 
is  based  on  the  maximum  range  of  the  radar  operations  R  ,  such  that 


f 


r 


c/2R 


max 


The  wavelength  is  then 


A  =  2v  / f  K 

max  r 


and  the  dwell  time  required  for  all  p  bursts  in  one  beam  position  is 


T 


d 


k  S/f  K 
P  r 


so  that  the  number  of  pulses  in  the  basic  burst  is 

N  =  k  e/Kp 
P 

A  practical  design  procedure  is  actually  quite  different  from  the  above 
procedure  because  the  penalty  for  a  large  number  of  pursts  p  is  not  factored 
into  the  optimization.  First,  a  large  value  of  p  implies  a  small  number  of 
pulses  N  per  burst.  Since  the  optimization  does  not  recognize  that  N  is  only 
an  integer,  a  practical  application  of  the  results  causes  one  to  round  up  to 


2-11 


the  next  higher  integer  value  of  N.  If  N  is  small  then  this  rounding  up  can 
represent  a  large  increase  in  dwell  time  beyond  what  was  used  in  the  optimiza¬ 
tion.  Second,  a  large  value  of  p  implies  the  radiation  of  a  great  deal  of 
radar  power  just  to  solve  the  blind-speed  problem.  That  is,  p  times  more 
power  is  required  than  is  required  when  only  one  burst  is  used.  It  is  there¬ 
fore  important  to  recognize  that  a  natural  bias  away  from  the  optimum  and 
toward  lower  p  is  preferable.  An  easy  means  of  recognizing  the  impact  of  lower 
than  optimum  p  is  therefore  required. 

Just  such  a  format  for  the  graphical  solution  has  been  developed,  as  illus¬ 
trated  in  Figures2-7  and  2-8.  In  this  format  each  new  value  of  p  represents  a 
separate  graph,  and  a  family  of  curves  representing  a  variety  of  values  of  R 
is  shown  on  each  graph.  It  is  then  an  easy  matter  to  inspect  the  curve  for 
the  proper  g  and  to  see  not  only  the  optimum  (maximum  K)  but  also  how  much  K 
must  be  sacrificed  in  order  to  arrive  at  a  lower  p.  Since  K  is  proportional 
to  carrier  frequency  this  tradeoff  is  easily  interpreted  in  its  impact  on 
the  radar  system. 

While  the  above  graphical  solutions  and  formulae  are  sufficient  for 

sizing  the  radar  system,  they  do  not  define  the  actual  schedule  of  PRFs  pre- 

* 

cisely.  It  has  been  found,  however,  that  Jolmson's  tables  can  be  used  along 
with  Table  2-1  to  produce  the  actual  PRFs  to  a  reasonable  accuracy. 

When  many  different  parameters  and  designs  must  be  considered,  however, 
a  more  automated  approach  is  preferable.  We  therefore  describe  next  the 
computer  program  that  has  been  developed  for  automatically  solving  the  blind- 
speed  problem. 

2.2  AUTOMATED  DESIGN  PROCEDURES 
2 .2.1 _ Approach 

The  need  for  a  computer  implementation  of  the  waveform  optimization  i6 
clear  when  one  considers  that  a  wide  variety  of  environments  must  be  investigated. 
A  main  objective  of  this  project  is  to  determine  whether  a  radar  sensor  that 


*S.  M.  Johnson,  Multiple  Repetition  Frequency  Radar  Coverage  (Maximum 
No-Gap  Coverage) ,  Rand  Corp.  Report  R-969-PR.  August  1972. 
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Table  2-1.  Values  of  x  Corresponding  to  (K  »  p  ) 
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adapts  to  its  environment  has  a  major  advantage  over  a  nonadaptive  radar. 

This  adaptive  radar  must  therefore  be  able  to  reoptimize  its  waveforms  to 
environmental  changes.  Accordingly  either  the  radar  must  perform  a  reopti¬ 
mization  in  real  time  or  a  large  number  of  waveform  sets  must  be  precomputed 
and  stored  in  the  computer.  In  either  case  the  generation  of  a  large  munber 
of  waveform  designs  is  implied. 

In  automating  the  desigr  procedures  it  is  not  enough  to  merely  automate 
the  previously  described  graphical  procedures.  Instead  one  must  consider 
that  two  range  zones  must  be  covered  by  the  same  radar.  Thus,  when  a  minimum 
wavelength  that  is  achievable  in  one  zone  is  lower  than  the  minimum  wavelength 
that  is  achievable  in  the  other  zone,  the  former  minimum  cannot  be  used.  In¬ 
stead  the  longer  wavelength  must  be  considered  as  a  constraint  in  determining 
the  waveforms  for  the  former  zone. 

The  automated  design  procedures  must  therefore  be  separated  into  two  ap¬ 
proaches.  In  practice  one  may  first  determine  the  minimum  allowable  wavelength 
for  the  near-range  zone  and  then  the  minimum  for  the  far-range  zone.  The 
actual  wavelength  will  have  to  be  the  greater  of  these  two  minima.  Then  the 
waveform  for  the  zone  that  results  in  the  smaller  wavelength  will  have  to  be 
redesigned  with  the  greater  wavelength  as  a  constraint.  The  overall  result 
will  be  a  two-zone  design  that  provides  both  short  search  times  and  small  wave¬ 
lengths. 

The  first  design  approach  is  therefore  to  find  the  solution  that  provides 
the  minimum  wavelength  for  given  environmental  parameters  B  and  p.  We  shall 
hereafter  call  this  approach  "Case  1."  The  second  design  approach  is  to  find 
the  solution  that  provides  the  minimum  required  search  dwell  time  for  given 
wavelength  (and  hence  improvement  factor  K)  and  parameter  p.  We  shall  call 
this  approach  "Case  2."  The  computer  solution  for  each  approach  is  detailed 
next . 

2.2.2  Automated  Design  for  Fixed  Wavelength 

Formulation  of  "Case  2" 

We  shall  detail  "Case  2"  first  because  it  will  be  easier  to  appreciate 

the  details  of  "Case  1"  once  the  solutions  for  "Case  2"  are  at  hand.  "Case  2" 

is  the  design  approach  in  which  the  system  wavelength  has  already  been  selected 

and  now  we  want  to  find  the  radar  waveform  that  minimizes  the  required  search 

dwell  time  for  a  given  range  rate  extent  AR  ^  for  the  clutter  and  maximum  speed 

v  for  the  target, 
max 
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Expressed  mathematically,  we  want  to  minimize  Tg  for  a  fixed  X  and  p.  Re¬ 
call,  however,  that  wavelength  is  related  to  the  blind-speed  improvement  factor  K 
as 


A  =  2  v  /f  K 
max  r 


where  f  is  that  maximum  allowable  pulse  repetition  frequency  (PRF) 


f 

r 


=  c/2R 


max 


where  R  is  the  maximum  range  for  the  zone  under  design.  Thus,  since  A, 
max 

v  ,  and  f  are  fixed,  K  is  also  fixed.  Recall  that  the  parameter  6  is: 
max  r 


R  = 


k  v 
a  max 


T 

s 


\  p  a  s  / 

which  is  proportional  to  the  total  dwell  time  T^ .  Usually  the  parameters  within 
the  parenthesis  will  be  prespecified  so  that  they  may  be  considered  fixed  for  the 
purposes  of  the  automated  design.  The  entire  design  objective  for  "Case  2"  can 
therefore  be  stated  parametrically  as:  "minimize  B  for  a  given  K  and  p." 


Optimization  Approach 

The  optimum  solution  is  subject  to  two  design  equations  (as  described  in 
Howard,  et  al,  op  cit) : 


K  = 


xP  -  2x 


P-1  + 


K  < 


x  -  1 

26 


(3-1) 

(3-2) 


-  x(2p3  +  p) 

where  x  is  a  normalized  PRF  (maximum  allowed  PRF  divided  by  Doppler  spread  of 
the  clutter)  given  by 

x  =  Af  / 2AR  „ 
r  cX 


(Recall  also  that  K  and  p  must  be  integers.)  The  constraint  given  by  (2)  can 
be  written  also  as 


o  >  pKx  (K  ,  p) _ 

-2(1-  Kpx (K , p) ] 


(3-3) 
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where  we  have  explicitly  shown  the  dependence  of  x  on  K  and  p.  &  is  therefore 
minimum  when  equality  holds  in  (3)  and  the  right-hand  side  of  (3)  is  minimized. 
Since  K  and  p  are  prespecified,  however,  the  only  parameter  that  can  be  adjusted 
to  minimize  &  is  the  number  of  bursts  p.  We  will  therefore  proceed  to  describe 
the  optimization  for  the  range  of  parameters  of  interest.  (From  experience 
this  range  is  p  from  1  to  12  and  K  from  1  to  60.) 

If  K  is  unity  then  we  should  of  course  use  p  as  unity.  If  K  is  an  inte¬ 
ger  larger  than  one  then,  due  to  the  nature  of  x(K,p),  the  right-hand  side  of 
(3)  will  have  a  single  positive  minimum  with  respect  to  p  only  at  some  p  >  2. 
Thus  we  may  implement  a  computer  search  for  a  minimum  by  starting  at  p=2  and 
evaluating  (3),  then  increasing  p  one  at  a  time  and  reevaluating  (3)  until 
the  value  of  (3)  starts  increasing.  Then  the  minimum  of  (3)  will  be  recog¬ 
nized.  This  "brute  force"  approach  is  reasonable  because  there  are  at  most 
only  eleven  values  of  p  to  consider.  (Recall  that  the  transmission  efficiency 
is  only  1/p  when  p  bursts  are  required  for  blind-speed  resolution.  The  smaller 
values  of  p  are  therefore  preferred,  and  p  >  12  is  unacceptable.) 

In  evaluating  (3),  however,  we  need  not  compute  the  function  x(K,p)  anew 
eleven  times!  Since  only  59  values  of  K  and  11  values  of  p  are  ever  of  interest 
for  the  subject  radar,  we  may  precompute  the  corresponding  649  (59  times  11) 
possible  values  of  x(K,p)  and  store  them  in  a  table.  Then  the  evaluation  of  (3) 
is  reduced  to  a  small  amount  of  arithmetic  and  a  serial  table  reading.  Such  a 
table  has  been  generated  off-line  by  using  a  recursive  algorithm: 

xi+i=xi+AXi  1=1-2.  ... 

where 

x  P  -  2x  P_1  -  Kx  +  K  +  1 

Ax  =  —  ‘ - - 

pxjLP  -  2(p-l)xip  -  K 

and 

X,  =  |j/(p-l)  +  1/K1/<P-1> 

-12 

We  iterate  until  Ax,  <  10  and  use  the  resultant  value  at  x.  , 

l  i+1 

for  x(K,p). 


as  our  value 
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Sub_tl 

A  few  details  in  the  solutions  must  be  handled  as  follows.  For  smaller 
values  of  p  the  right-hand  side  of  (3)  can  be  negative  depending  upon  the  values 
of  K  and  p.  In  those  cases  there  is  no  value  of  6  which  satisfies  the  condi¬ 
tions.  However,  for  a  fixed  K  and  p  as  soon  as  a  non-negative  value  of  8  ex¬ 
ists  it  will  be  non-negative  for  any  larger  value  of  p.  Thus  we  can  determine 
the  minimum  p  to  be  used  to  start  the  search  for  the  minimum  &• 

One  can  show  that  if  p  >  1  or  K  >  l/2p  then  no  solution  exists.  On  the 

computer  output  we  shall  denote  this  condition  by  assigning  a  value  of  p=0  and 
3=0.  If  the  optimum  value  of  p  is  larger  than  12  we  shall  denote  this  condi¬ 
tion  by  assigning  p=13,  8=10^. 

If  K  is  1  then  p  is  1.  Under  these  circumstances  x  is  undetermined  in  (1) 

and  can  have  a  value  between  zero  and  one.  We  use  the  value  x=l  in  the  right- 

hand  side  in  this  case  and  thus  have  the  computer  output  p=l ,  8=l/[ 2 (1 -p)  ] . 

Subroutine  OPB 

Using  the  above  procedure  a  FORTRAN  subroutine  called  OPB  (K,p,Bm>p)  has 
been  generated.  For  a  given  K  and  p  it  generates  the  minimum  B(  =  8  )  and  the 
corresponding  p.  The  flow  diagram  for  this  subroutine  is  shown  in  Figure  2-9. 

2.2.3  Automa ted  Design  for  Fixed  Dwell  Time 
Optimization  Approach 

We  now  consider  "Case  1"  where  dwell  time  is  fixed  and  wavelength  must  be 
minimized.  In  this  case  B  and  p  are  given  and  we  want  the  maximum  integer  K 
along  with  its  associated  p  that  satisfy  the  design  equations  (1)  and  (2).  Since 
K  is  an  integer  it  is  possible  that  more  than  one  p  will  provide  the  same  inte¬ 
ger  value  for  K.  In  that  case  we  choose  the  smallest  of  the  alternative  values 
of  p  because  small  p  corresponds  to  high  transmission  efficiency. 

Another  way  of  stating  the  problem  is  that  we  wish  to  find  the  integer  K 
which,  for  the  given  p,  yields  a  B  (and  associated  p)  as  close  as  possible  to 
the  given  B-  (The  reason  the  problem  can  be  stated  in  this  form  is  the  fol¬ 
lowing.  A  given  K  and  p  will  have  a  minimum  8  and  corresponding  p.  Then 

these  values  of  8  and  p  will  give  a  non-integer  value  for  the  maximum  K.  This 
is  why  we  use  the  largest  minimum  8  which  does  not  exceed  the  actual  B  and 
yet  gives  an  integer  K.) 


K  and  p  are  inputs  to  the  subroutine.  The  minimum  3  (denoted  3  )  and 

m 

its  associated  p  are  outputs. 


Figure  2-9.  Subroutine  OPB  (K,p,B  »p) 
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Subroutine  OPK 

A  subroutine  to  compute  the  maximum  K  has  been  developed,  and  is  called 

OPK  (8,p,K,p,8  ).  (See  Figure  2-10.)  For  inputs  B  and  p  it  generates  the 
m 

maximum  integer  K  and  corresponding  integer  p.  It  also  generates  B^,  which 
is  the  value  of  B  needed  to  obtain  the  integer  value  of  K.  That  is,  ^ 

B  yields  as  high  an  integer  value  for  K  as  B>  and  at  the  same  time  the  p  cor- 

ID 

responding  to  f  is  the  smallest  value  which  will  yield  this  integer  value 

m 

for  K. 

The  subroutine  OPK  (B,P,K,p,B  )  uses  the  subroutine  OPK  (K,p,B  ,p)  in  the 

m  m 

following  way.  It  can  be  shown  from  (3)  and  the  nature  of  x(K,p)  that  B  in- 

ID 

creases  as  K  increases  for  a  fixed  p.  We  can  therefore  determine,  by  in¬ 
creasing  K  one  by  one,  which  K  yields  Bffi  closest  to  B  but  no  larger  than  B. 

The  subroutine  OPB  also  determines  the  corresponding  p. 

This  then  gives  us  the  maximum  value  of  K  and  the  corresponding  p.  Since 

B  is  a  minimum  and  in  general  less  than  B»  the  form  of  (3)  then  implies  that 
m 

lowering  p  will  increase  B  above  Bffi  for  K  fixed  at  its  maximum  integer  value. 

Thus  for  the  fixed  K  and  p  we  test  lower  values  of  p  to  find  the  smallest  p 
for  which  the  8  is  no  larger  than  the  B  given  in  OPK. 

Subroutine  SOPK 

This  is  done  using  the  equality  sign  in  (3)  in  a  subroutine  called  SOPK 
(B.K.p.p.B  )•  (See  Figure  2-11.)  The  inputs  to  SOPK  are  the  given  B  and  p 

TT1 

along  with  the  maximum  integer  K  and  corresponding  p  and  B  found  from  OPK. 

ID 

The  outputs  from  SOPK  are  8,K,p  and  the  smallest  p  and  associated  B  (still 
called  B  )  which  yield  the  maximum  integer  K  and  yet  has  B  <  B- 

ID  ID 

Since  the  value  of  K  that  we  are  seeking  lies  in  the  range  of  1  to  60,  it 
is  more  efficient  to  perform  a  binary  search  for  K  than  to  search  sequentially 
starting  at  K=1  and  proceeding  upwards  using  OPB  (K,p,B  ,p).  From  (1)  and  (2) 

ID 

we  can  show  that 

K  <  l/2p 

so  that  the  largest  value  to  use  for  K  in  the  binary  search  is  bounded.  The 
lowest  value  of  K  we  use  in  the  binary  search  is  2.  If  the  search  indicates  the 
resultant  K  is  less  than  2  we  take  K=l.  If  the  search  indicates  the  resultant  K 
is  larger  than  60  wc  use  K~6Q.  Again  p  is  assumed  no  larger  than  12  so  that  a 
p-12  will  result  if  in  reality  p  >  12.  In  this  case  the  value  of  K  obtained  from 
OPK  may  be  smaller  than  the  maximum  K. 
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8  and  p  are  inputs  to  subroutines.  The  maximum  integer  K  and  cor 

responding  p  are  outputs.  Also  generated  is  the  actual  6  <6 

tn  * 

needed  to  achieve  this  K. 

Figure  2-10.  Subroutine  OPK  ) 

m 
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If  p=l  or  2  then  return 


For  pn=p-l,  p-2,  ...  compute 
PnKx(K,pn) 

Bn  “  2[ l-Kpx(K,p  )] 
n 

until  8  >  6  or  8  is  negative, 

n  n 

Then  the  previous  p^  along  with 

the  previous  8  are  the  minimum 
n 

p  and  corresponding  6  <  8  which 

IS 

achieves  K.  For  pn=p-l  the  pre¬ 
vious  p  and  8  are  the  input 
n  n 

values  of  p  and  8  •  Return 


Figure  2-11.  Subroutine  SOPK  (8,K,p,p,8  ) 

IB 
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As  with  OPB  we  have  no  solution  if  p  >  1.  We  denote  no  solution  by  pro* 

ducing  K“0,  p“0  and  8  •  6  in  OPK.  Also,  from  the  form  of  (3)  and  the  fact  that 

m 

x  <  1  (even  if  K*l)  ,  we  can  show  that  no  solution  exists  when  8  <  1/(2 ( 1  —  p) J - 

If  the  maximum  allowable  K  is  1  then  the  corresponding  p  is  1.  We  then 
set  8  «  1 / [ 2 ( 1 — p) ]  because  any  8  larger  than  this  will  be  large  enough  to  pro- 

ID 

duce  K**l ,  p«l . 

2.2.4  Impact  of  Lower  than  Optimum  p 

OPK  (0.P.K.P.8  )  yields  the  maximum  integer  K  and  corresponding  p  for 

a  given  8  and  p.  Also  OPB  (K,p,8  .p)  yields  the  minimum  8=8  and  corresponding 

ID  ID 

p  for  a  given  K  and  p.  In  these  two  optimization  problems  it  is  of  interest  to 
know  how  much  is  lost  in  the  optimizations  if  smaller  values  of  p  than  the  opti¬ 
mum  are  used.  Due  to  average  energy  constraints  and/or  processor  complexity 
there  are  designs  where  it  is  advantageous  to  lower  the  value  of  p  despite  a 
sacrifice  from  optimality. 

Thus  given  3,p  and  p  (less  than  the  optimum  p)  we  wish  to  solve  for  K  after 
we  have  used  OPK.  Given  K,  p  and  p  (again  less  than  the  optimum  p)  we  wish  to 
solve  for  8- 

The  former  problem  can  be  most  easily  solved  by  using  (2).  Substituting 
8,  p  and  the  p  of  concern  into  the  right-hand  side  of  (2)  we  put  decreasing 
integer  values  of  K  into  the  right-hand  side  (beginning  with  the  optimum  K) 
until  we  reach  a  K  where  the  inequality  holds.  The  resulting  K  is  the 
achievable  K. 

This  suboptimum  K  is  found  by  using  the  defined  function  called 

FK(8«P»P  »K  ,p).  The  flow  chart  for  this  function  is  shown  in  Figure  2-12. 
o  o 

The  inputs  are  the  values  of  8  and  p.  •  Also  required  are  the  resultant  optimum 
values  of  p  and  K  (called  p^  and  K^)  and  the  value  of  p  for  which  we  wish  to 
calculate  K  (called  FK) . 

For  the  latter  problem  we  can  simply  calculate  8  using  the  equality  sign 
in  (3)  (i.e.,  K,  p  and  p  are  given).  We  do  thi6  using  the  FORTRAN  function 
FB(K,p,p)  shown  in  Figure  2-13. 

If  p  is  made  too  small  we  may  have  a  situation  where  8  does  not  exist.  In 
this  case  the  right-hand  side  of  (3)  is  negative.  If  this  situation  occurs  we 
output  FB=0  in  Figure  2-13. 

In  Appendix  A  we  have  listed  the  code  for  the  subroutines  OPB  (K,p,B  ,p)  , 

ID 

OPK  (8,p>K,p,B  )  and  SOPK  (8.K,p,p,8  )  as  well  as  the  functions  FK  (B,p,p  ,K  ,p) 
m  m  o  o 

and  FB  (K,p,p). 
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2.3  SPECIFIC  BASELINE  DESIGNS 

In  order  to  illustrate  the  use  of  the  new  automated  design  procedures  we 
use  here  the  baseline  parameters  from  our  first  report  (Howard,  et  al. ,  op  cii)  . 
Two  distinct  designs  based  on  two  different  antenna  concepts  were  developed. 

The  first  concept  used  a  set  of  four  rotating  antennas  with  electronic  scan  in 
elevation.  The  second  concept,  which  is  clearly  more  costly  but  also  more  capable 
in  tracking,  used  a  set  of  four  stationary  antennas  with  electronic  scan  in  both 
elevation  and  azimuth. 

2.3.1  Design  for  Rotating  1-D  Phased  Array 
Background 

For  the  mechanically  rotating  antenna6  we  are  greatly  constrained  in  the 
use  of  the  antennas  for  various  tasks.  Each  of  the  four  antennas  must  be  used 
partly  for  tracking  in  order  to  achieve  the  required  2  Hz  tracking  rate.  (The 
antenna  rotation  rate  is  30  rpm.)  Roughly  half  of  the  time  load  for  each  an¬ 
tenna  face  is  devoted  to  tracking.  For  one  of  the  antenna  faces  the  other  half 
of  the  load  is  reserved  for  long-range  search  (search  in  the  far-range  zone) 
in  which  the  frame  time  is  allowed  to  be  as  many  as  5  or  so  antenna  revolutions. 
For  the  other  three  faces  the  other  half  of  the  load  is  reserved  for  short- 
range  search  (in  the  near-range  zone)  in  which  the  frame  time  is  only  2  seconds 
so  that  all  bursts  required  for  the  blind-speed  resolution  must  be  implemented 
by  the  three  faces  in  just  one  revolution  of  the  antenna. 

Far-range  Zone 

Let  us  consider  the  far-range  zone  first.  For  search  in  this  zone  we 
use  just  one  antenna  and  devote  only  half  of  its  time  to  search.  Thus,  in 
every  antenna  revolution  we  have  available  only  1  second  for  search.  However, 
we  allow  the  frame  time  to  be  2  p  seconds  where  p  is  the  number  of  bursts  re¬ 
quired  for  blind-speed  resolution,  as  long  as  p  does  not  exceed  about  3  or  so. 

Thus  the  total  search  time  available  for  this  mode  is  p  seconds  for  a  search 
over  2n  in  azimuth.  We  also  have  no  ground  clutter  in  the  long-range  search, 

:.<>  that  the  clutter  .spiend  is  just  the;  range-rate  spread  of  wind-driven  rain, 
or  5  mps.  The  parameters  for  the  long-range  search  are  therefore 

k  =0.85 
a 

k  =1.0 
P 

v  =  1000  mps 
max  r 


i 


i 
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w  =  3  m 
a 

T  =  p  sec 
s 

e  =  2v 

s 

A*cJt  =  5  mps 

so  Chat  the  inputs  to  the  computer  design  programs  are 

B  =  45.1  p 

p  =  0.005 

The  problem  could  now  be  stated  that  we  wish  to  find  the  optimum  K,  p  for 
a  given  p,  B  but  where  B  is  dependent  on  the  solution  for  p.  Clearly,  if  B 
were  fixed,  we  could  use  the  programs  to  find  the  optimum  K,  p.  However,  since 
B  is  proportional  to  p  a  different  procedure  is  required.  Recall  that  only  a 
few  values  of  the  integer  p  are  of  interest.  We  could  therefore  approach  the 
problem  by  finding  K  for  a  given  p,  B»  and  p,  first  for  p=l  then  repeating  the 
procedure  for  all  integer  p's  of  interest.  We  could  then  choose  the  K,  p  pair 
that  is  optimum  or  at  least  acceptable.  This  procedure  is  performed  by  re¬ 
peatedly  using  subroutine  OPK  and  function  FK.  For  a  preselected  value  of  p  we 
compute  B  (as  45.1  p  in  the  example)  and  the  use  B,  P  as  inputs  to  subroutine 
OPK.  The  output  of  the  subroutine  is  the  optimum  pair  K^,  p^  for  the  given 

fixed  values  for  B»  P-  Thi6  optimum  pair  is  now  used,  along  with  the  pre¬ 

selected  values  for  p,  B,  and  p,  as  the  input  to  the  function  FK.  The  value 
of  the  function  will  then  be  the  value  of  K  for  the  given  values  of  p,  B  and 
P. 

When  the  above  procedure  is  followed  for  the  parameters  of  the  example. 
Table  2-2  and  Figure  2-14  result.  Evidently  the  achievable  improvement  fac¬ 
tor  K  increases  without  limit  as  the  number  of  bursts  p  is  increased!  In 

other  words,  by  not  constraining  the  search  time  and  instead  allowing  it  to 

grow  with  the  number  of  bursts  used  for  blind-speed  resolution  we  find  the 
very  high  improvement  factors  (and  hence  short  wavelengths)  are  achievable 
without  limit. 

Evidently  this  conclusion  is  unrealistic.  We  had  previously  stated  that 
we  allow  the  search  frame  time  to  increase  with  p.  A  high  p  therefore  implies 
not  a  need  for  more  transmitted  power  but  a  need  for  a  long  frame  time.  The 


Figure  2-14.  Example  of  Optimization  for  Long-Range  Search 
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frame  time  is  2  p  seconds  for  the  example.  An  increased  frame  time  is  there¬ 
fore  acceptable  for  the  far-range  zone  as  long  as  p  is  not  greater  than  5.  By 
contrast,  frame  times  longer  than  10  seconds  would  probably  not  be  acceptable 
because  of  the  possible  extensive  maneuvers  of  the  targets  within  such  a  long 
frame.  Thus  p  <  5  is  a  realistic  constraint  for  the  example. 

If  we  choose  p=4  we  find  from  Figure  2-14  that  the  achievable  improvement 
factor  is  19.  This  solution  results  in  an  8-second  frame  time,  as  desired. 
Since  the  maximum  range  for  the  far-range  zone  is  200  km,  the  maximum  allowable 
PRF  for  the  bursts  is 

f  -  c/2  R  -  750  Hz 
r  max 

Accordingly  the  minimum  achievable  wavelength  is 

X  ■=  2  v  /f  K  =  14.0  cm 
max  r 


Near- range  Zone 

In  the  near-range  zone  we  are  even  more  constrained  because  we  must  per¬ 
form  the  search  within  a  single  revolution  of  the  antennas.  Recall  that  one 
of  the  four  antennas  has  been  reserved  for  the  search  of  the  far-range  zone. 
The  remaining  three  antennas  are  therefore  potentially  available  for  the 
search  of  the  near- range  zone.  However,  half  of  the  available  time  resources 
is  reserved  for  tracking.  Therefore  each  antenna  has  available  for  near¬ 
range  search  a  1-second  dwell  time  out  of  the  specified  2-second  near-range 
frame  time.  Since  three  antennal  are  operating  simultaneously  fo’-  performing 
the  search,  it  is  possible  to  transmit  separate  bursts  on  separate  antennas 
at  the  same  time.  On  the  other  hand,  since  coherency  must  be  maintained 
within  each  burst  so  that  Doppler  filtering  is  provided,  the  dwell  time  per 
burst  is  limited  to  2  second  for  1,  2,  or  3  bursts  (p=l,2,  or  3),  0.5  seconds 
for  4,5,  or  6  bursts,  0.33  seconds  for  7,8,  or  9  bursts,  etc.  Evidently  the 
available  time  is  not  efficiently  used  when  p  is  not  a  multiple  of  3.  That  is, 
the  total  dwell  available  is  really  3  seconds  because  there  are  3  antennas. 

But  only  2  seconds  are  used  when  p=2,  2.5  seconds  are  used  when  p=5 ,  etc., 
because  of  the  constraint  on  dwell  time  per  burst.  Thus,  for  the  mechani¬ 
cally  rotating  antenna  we  are  biased  toward  the  solutions  for  p=3,  6,  and  9. 
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Keeping  in  mind  this  constraint  on  p  we  may  now  attempt  to  achieve  a 
satisfactory  design  as  follows.  As  far  as  the  mathematics  and  computer  pro¬ 
grams  are  concerned  we  can  now  develop  the  design  by  considering  one  antenna 
searching  a  2?f/3  azimuth  sector  rather  than  three  parallel  antennas  searching 
a  2tt  sector.  The  relevant  design  parameters  are  therefore: 

k  =0.85 
a 

k  =1.0 
P 

v  =  1000  mps 
max 

w  =  3  m 
a 

T  =1.0 
s 

6  =  2tt/3 

s 

=  20  mps  (including  ground  clutter) 

so  that  the  inputs  to  the  computer  subroutine  are 

B  =  135.3  (equivalent  to  a  =  270.6) 
p  =  0.02 


The  optimum  solution  that  results  from  an  application  of  subroutine  OPK 
is  K=10,  p=5.  However,  a  nearly  optimum  solution  exists  at  K=9,  p=4.  Further¬ 
more,  in  light  of  the  3-antenna  constraint  a  better  solution  is  actually 
K=8,  p=3.  The  maximum  range  in  the  near-range  zone  is  80  km  so  that  the  PRF 
is  nominally 

f  =  c/2R  =  2875 
r  max 

Accordingly,  the  minimum  allowable  wavelength  is 

A  =  2  v  /f  K  =  13.3  cm 
max  r 

(Since  we  have  not  selected  the  mathematical  optimum  the  actual  dwell  time 
needed  for  searc..  may  be  less  than  the  originally  allowed  time.  In  order  to 
find  the  actual  time  now  needed  we  may  recompute  the  parameter  B  for  the  se¬ 
lected  values  of  K  and  p.  When  we  do  this  we  find  that  B  is  98  rather  than 
135.3.  The  time  needed  for  the  near-range  search  is  therefore  only  73%  of 
the  originally  allocated  time  of  1  second  per  antenna  per  search  frame.) 
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The  13.3  cm  wavelength  o£  this  solution  is  slightly  shorter  than  the 
minimum  allowable  wavelength  for  the  far-range  zone.  It  is  therefore  evi¬ 
dent  that  the  far-range  zone  determines  the  operating  frequency  in  this  ex¬ 
ample.  Thus,  we  must  now  reoptimize  the  search  in  the  near-range  zone  by 
relaxing  the  wavelength  to  A  =  14.0  cm  and  then  solving  for  the  actual  al¬ 
lowable  search  time.  We  expect  to  find  that  a  somewhat  lower  6  is  optimum 
for  p=3 ,  thereby  allowing  a  reduction  in  the  time  that  must  be  devoted  to 
the  short-range  search.  This  reduction  would  make  available  more  time  for 
tracking,  thereby  improving  the  design  beyond  what  would  result  if  we  merely 
retained  the  bursts  determined  for  13.3  cm  wavelength. 

When  we  attempt  the  reoptimization,  however,  we  find  that  a  relaxation 
of  wavelength  from  13.3  cm  to  14  cm  is  not  enough  to  materially  alter  the  solu¬ 
tion.  In  fact,  we  recall  that  K  is  related  to  A  through  the  parameters  f  and 

v  so  that,  for  A=14  cm  we  have 
max 


K 


2v 

max 


/Af 


r 


7.6 


Thus,  since  A  has  not  been  increased  sufficiently  to  change  K  by  a  full  integer, 
we  must  retain  the  preceding  solution  of  K=8,  p=3.  However,  we  may  improve 
the  design  by  adjusting  the  maximum  PRF  downwards  such  that  K  is  exactly  equal 
to  8.  That  is,  we  select  f  =  1786  Hz  (recalling  that  v  is  1000  mps ,  A  is 

IT  IT13X 

14  cm).  The  range  ambiguities  are  now  at  a  farther  range  than  previously  so 
that  the  maximum  range  of  the  near-range  zone  can  be  increased  to 

R  =  c/2f  =  84  km 
max  r 

This  increase  in  maximum  range  is  the  only  practical  benefit  from  the  slightly 
lowered  wavelength. 

Swenaiy  of  the  Deuicjn 

The  baseline  design  for  the  mechanically  rotating  set  of  antennas  can  now 
he  summarized  from  the  preceding  results  along  with  .lohnson's  tables.  The  far- 
range  zone  is  searched  with  one  antenna  in  four  successive  revolutions  of  the 
antenna.  A  different  PRF  is  used  in  each  revolution  in  order  to  resolve  all 


( 


2-33 


blind  speeds  within  the  four  revolutions  (or  8  seconds).  The  PRFs  for  the 
solution  <K=19,  p=4,  x=3.23)  are,  using  an  automatic  algorithm  for  computing 

A 

PRFs  (to  be  described  later) 


f  =  750  Hz 
rl 

f  „  =  722.75  Hz 
r2 

f  =  668.26  Hz 
r3 

f  ,  =522.94  Hz 
r4 

(This  set  of  PRFs  leads  to  just  one  blind  region  at  a  Doppler  frequency  of  63f  Hz 
and  5.25  Hz  wide.)  The  number  of  pulses  at  each  PRF  are 


so  that  the  total  dwell  time  for  all  4  bursts  in  one  beam  position  is  28.7  milli¬ 
seconds  (rather  than  the  theoretical  T.  =  25.3  milliseconds).  Thus  56.5%  of  the 

a 

time  load  for  one  antenna  is  devoted  to  the  long-range  search. 

The  near-range  zone  is  searched  with  three  antennas  in  one  revolution  of 
the  antenna.  A  differeit  PRF  is  used  for  each  successive  look  at  the  search 
area.  The  PRFs  for  the  solution  (K=8,  p=3,  x=3.54)  are,  using  the  automatic 
algorithm  once  again 


rl 


1786  Hz 


f  2  =  1625.79  Hz 

f  .  =  1345.42  Hz 
r3 

(This  solution  produces  a  single  63.75  Hz  wide  blind  zone 
The  number  of  pulses  at  each  PRF  are 
n  =  12 
n2  =  11 

n3  =  9 

so  that  the  total  dwell  time  for  all  3  bursts  in  one  beam 

(rather  than  the  theoretical  T,  =  18.9  msec).  Thus  53.4% 

d 

each  antenna  is  devoted  to  the  near-range  search  zone. 


centered  at  1566  Hz.) 


position  is  20.2  msec 
of  the  time  load  for 


At  this  point  we  could  have  used  Johnson's  tables,  as  described  in  our 
earlier  report.  However,  when  we  developed  techniques  for  determining  the  set 
of  PRFs  in  real  time,  we  found  a  better  design  approach,  as  detailed  in 
Section  4. 
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2.3 .2  Design  for  Stationary  2-D  Phased  Array 
Background 

We  are  much  less  constrained  with  a  2-D  phased  array  in  that  the  same 
antenna  provides  all  bursts  so  that  we  have  complete  flexibility  in  the  num¬ 
ber  of  bursts  allowed.  Thus  we  need  not  define  allowable  search  times  consis¬ 
tent  with  the  periodic  reappearance  of  rotating  antennas.  If  we  assume  a  total 
frame  time  in  the  far- range  zone  of  8  sec,  then  we  may  perform  the  design  for 
a  total  dwell  of  1.5  seconds  for  a  search  over  a  7r/2  sector.  (The  other 
three  stationary  antennas  are  designed  identically.)  Similarly,  in  the  short- 
range  zone  we  no  longer  find  that  p=3  and  p=6  are  the  only  easily  implemented 
solutions.  Thus  such  solutions  as  p=4  are  allowed. 

Far-range  Zone 

For  the  long-range  search  we  now  have  the  parameters: 


k  =  0.85 
a 

k  =  1.0 
P 

v  =  1000  mps 
max 

w  =  3m 
a 

T  B  1.5  sec  (for  8  sec  frame  time) 
s 

e  =  n/2 

s 

AR^  =  5  mps 

so  that  the  inputs  to  the  computer  programs  are 


6  =  270.6 
p  =  0.005 

The  computer  solutions  (or  graphical  solutions  from  Figure  2-5)  are  K=25,  p=5, 

x=2.75.  However,  we  also  find  little  sacrifice  in  K  by  choosing  p=4,  for 

which  K=23.  Moreover,  a  25%  increase  in  power  and  time  efficiency  results, 

and  the  longer  attendant  wavelength  makes  the  elevation  beam  sufficiently 

wide  that  a  major  gain  loss  from  cosecant-squared  beam  shaping  does  not  result. 

The  PRF  and  wavelength  corresponding  to  the  solution  are 

f  =  c/2R  =  750  Hz 
r  max 

X  =  2v  /f  K  =  11.6  cm 
max  r 

Kvidi-ntly  we  have  gained  a  great  deal  in  terms  of  reducing  X  merely  by  allowing 
more  overall  time  for  the  long-range  search  because  the  constraints  of  the  mech¬ 
anical  array  have  been  removed.  (Because  we  do  not  have  to  rely  on  revolving 


antennas  to  achieve  the  required  2  Hz  tracking  rate,  we  no  longer  have  to  re¬ 
strict  the  search  time  to  half  the  load  of  one  antenna,  or  1.0  seconds  per 
revolution  per  antenna.  Thus  the  increase  in  allowance  to  1.5  seconds  per  an¬ 
tenna  is  a  consequence  of  the  removal  of  the  mechanical  constraint.)  Overall, 
18.75%  of  the  time  load  for  the  antenna  is  devoted  to  long-range  search. 

Near  Range  Zone 

In  the  near-range  zone  we  may  work  backwards  from  the  design  wavelength. 
That  is,  for  the  80  km  maximum  range  in  the  near-range  zone  we  have  the  maxi¬ 


mum  PRF  of 


f  =1875  Hz 
r 


Thus  we  have  for  the  selected  11.6  cm  wavelength  a  value  for  improvement  facto 


K  =  2v  /f  X  =  9.2 
max  r 

Since  X  is  fixed  at  11.6  cm,  we  must  use  an  integer  value  of  K  larger 

than  9.2.  Otherwise  v  would  have  to  decrease  and/or  f  would  have  to  in- 

max  r 

crease.  Neither  situation  is  desirable.  On  the  other  hand  we  do  not  wish  to 
use  a  K  any  larger  than  necessary  as  the  larger  K,  the  larger  a  so  that  the 
scanning  time  would  increase  unnecessarily.  Thus  we  take 

K  =  10 

Since  10  is  larger  than  9.2  we  can  lower  the  PRF  and  thus  extend  the  maximum 
range  in  the  short-range  mode.  The  result  is 


=  1724  Hz 


The  corresponding  maximum  range  is 


R  =  — -  =  87  km 

max  2f 

r 

The  value  for  K  of  10,  coupled  with  the  value  for  p  of  0.02  in  the  near-range 
zone,  leads  to  an  optimum  solution  for  B  and  p.  The  result  is  8=115,  p=5 . 

The  amount  of  time  devoted  to  the  short-range  search  is  therefore  (for  one  an¬ 


t  enna) 


BkpWa6s  115  x  1  x  3  x  n/2 


k  v 
a  max 


.85  x  1000 


0.638  sec 
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out  of  the  available  2-second  frame  time.  Thus,  only  31.9%  of  the  total  time 
for  each  antenna  need  be  devoted  to  the  short-range  search. 


S ternary  of  the  Design 

In  summary  we  use  18.75%  of  the  time  for  long-range  search,  31.9%  of  the 
time  for  short-range  search,  leaving  49.35%  of  the  time  for  tracking.  More¬ 
over,  this  has  been  done  with  a  wavelength  of  just  11.6  cm.  This  contrasts 


with  the  mechanical  antenna  set  where  less  than  50%  of  the  time  remains  for 


tracking  and  the  wavelength  is  14.8  cm! 

We  may  now  determine  the  actual  schedule  of  bursts  to  be  used  for  blind- 
speed  resolution.  The  set  of  PHI'S  for  the  long-range  search  is  found  from 


our  new  algorithm  as 


f  =  750  Hz 
rl 

f  „  =  727.02  Hz 
r2 

f  ,  =  676.46  Hz 
r3 

f  .  =562.70  Hz 
r4 


where  a  single  blind  zone  of  34  Hz  width  remains.  The  number  of  pulses  at 
each  PRF  are 


n^  =  5 


From  these  results  we  find  a  total  dwell  time  for  all  4  bursts  of  34.7  msec 
(rather  than  the  theoretical  T^  =  31.4  msec).  In  this  case  20.7%  of  the  time 
load  for  one  antenna  is  devoted  to  the  long-range  search. 

The  set  of  PRFs  for  the  short-range  search  is 

f  ,  =  1724  Hz 
rl 

f  ,  =  1623.61  Hz 
rl 

f  ,  =  1510.68  Hz 
r3 

fr4  =  1313.05  Hz 

f  =  1016.60  Hz 
r5 

In  this  case  no  blind  zones  remain.  The  number  of  pulses  at  each  PRF  are 
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These  results  lead  to  a  total  dwell  time  for  all  5  bursts  of  14.6  msec  (rather 

than  the  theoretical  T,  =  13. 3  msec ,  Then  35  %  of  the  time  load  for  each  antenna 

d 

is  devoted  to  the  short-range  search. 
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3.  OTHER  DESIGN  APPROACHES 

Two  other  design  approaches  were  also  studied  in  a  further  search  for 
extended  performance.  The  first  approach  is  a  variation  of  the  baseline 
blind-speed  resolution  scheme  in  which  multiple  bursts  are  transmitted  simul¬ 
taneously  on  different  frequencies.  The  main  question  in  its  application  is 
whether  the  added  performance  warrants  the  added  equipment  required.  The 
second  approach  is  an  extension  of  the  baseline  scheme  to  the  resolution  of 
blind  ranges  as  well  as  blind  speeds.  The  main  question  here  is  whether  the 
added  performance  warrants  the  added  complexity  and  sidelobe  requirements. 

3.1  MULTIPLE  SIMULTANEOUS  FREQUENCIES 

3.1.1  Basic  Principles 

We  have  seen  that  the  successful  resolution  of  blind  speeds  consumes  a 
considerable  amount  of  extra  dwell  time.  Our  past  approach  has  always  been 
to  transmit  sequentially  the  p  bursts  that  are  required.  However,  if  we  can 
provide  isolated  parallel  transmitting  channels  at  different  frequencies  then 
we  can  conceivably  transmit  the  bursts  in  parallel  and  thereby  save  a  consider¬ 
able  amount  of  dwell  time.  Variations  of  this  approach  are  detailed  here. 

The  parallel  use  of  multiple  carrier  frequencies  in  place  of  the  sequentia 
use  of  a  single  carrier  frequency  is  acceptable  in  principle  as  long  as  one 
recognizes  that  it  is  range  rate,  not  Doppler  frequency,  that  must  be  made 
blind-zone  free.  Since  range  rate  is  related  to  Doppler  frequency  as 

R  =  Afd/2  , 

the  ambiguous  range  rate  is  related  to  PRF  as 


Thus  a  change  in  wavelength  is  equivalent  to  a  change  in  PRF  as  far  as  blind- 
speed  resolution  is  concerned.  Moreover,  any  new  wavelength  that  is  suffi¬ 
ciently  far  for  isolation  from  the  basic  wavelength  could  be  used  as  long  as 
the  PRF  were  adjusted  accordingly.  The  use  of  multiple  frequencies  therefore 
provides  an  easy  means  of  cutting  the  dwell  time  to  reasonable  levels. 

The  conceptually  most  simple  implementation  would  be  to  use  p  carrier 
frequencies  with  wavelengths  proportional  to  the  multiple  PRFs  that  were 
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formerly  required  for  blind-speed  resolution.  With  this  concept  the  s.une 
pulsing  and  PR F  would  be  used  for  all  carrier  fiequencies  simultaneously. 

The  approach  therefore  appears  simple  and  attractive.  However,  a  number  of 
shortcomings  become  clear  upon  further  examination.  First,  the  range  of  car¬ 
rier  frequencies  required  is  usually  too  large  for  most  radars.  For  example, 
the  PRFs  required  for  the  baseline  design  for  the  rotating  antenna  differ  by 
factors  up  to  1.5.  This  means  that  the  multi-frequency  approach  would  re¬ 
quire  a  very  wide  range  of  operating  frequencies,  such  as  2  GHz  to  3  GHz. 
Second,  the  tuning  of  the  multiple  frequencies  is  rather  critical  for  success¬ 
ful  elimination  of  blind  zones.  This  means  that  when  the  environment  changes 
such  that  a  new  set  of  PRFs  would  conventionally  be  required,  the  multiple  fre¬ 
quencies  must  be  precisely  retuned.  A  truly  adaptive  radar  would  therefore  bt 
considerably  more  difficult  to  implement  with  multiple  frequencies  than  with 
multiple  PRFs. 

The  more  reasonable  approach  from  the  standpoint  of  equipment  is  the  use 
of  separate  carrier  frequencies  for  isolation  only ,  and  the  use  of  PRF  ad¬ 
justments  for  the  actual  blind-speed  resolution.  With  this  approach  there 
are  actually  several  levels  of  sophistication.  For  example,  one  could  use 
a  two-channel  system  in  which  two  PRFs  are  sequentially  and  independently 
transmitted  in  each  channel  such  that  the  blind-speed  resolving  capabilities 
of  a  four-PRF  system  are  achieved.  Or  one  could  use  a  four-channel  system 
with  just  one  PRF  for  each  channel  (albeit  each  PRF  is  different).  We  there¬ 
fore  can  achieve  different  levels  of  dwell  time  reduction  depending  on  the 
different  number  of  parallel  channels  available. 

Designs  for  a  multi-frequency  radar  can  be  developed  with  the  aid  of 
the  automatic  design  procedures  as  follows.  First  it  must  be  assumed  that 
the  multiple  frequencies  do  not  differ  in  wavelength  by  more  than  about  10%. 
With  this  assumption  the  change  in  Doppler  resolution  with  change  in  wave¬ 
length  can  be  ignored.  Second,  the  basic  (highest)  PRF  will  be  used  at  the 
lowest  carrier  frequency  (longest  wavelength).  Denoting  this  PRF  and  wave¬ 
length  as  i TO  and  ,  we  first  design  the  system  for  this  wavelength  alone. 

However,  we  replace  the  total  dwell  l  ime  available  for  search,  T  ,  by  the 

s 

quantity  qT^  where  q  is  the  number  of  parallel  channels  available.  (Note 
that  q  <  p  is  required,  and  furthermore  it  is  assumed  that  p/q  is  an  integer 
so  that  all  channels  are  used  to  their  full  capacity.)  The  set  of  PRFs 


3-3 


that  results  from  a  blind-speed  optimization  is  then  modified  by  the  wave¬ 
length  of  the  available  channels.  In  general  the  lowest  PRF  is  used  with 
the  shortest  wavelength  so  that  the  adjusted  PRF  is  raised  toward  the  basic 
PRF.  In  other  words,  the  multiple  wavelengths  available  will  have  a  nor¬ 
malizing  effect  on  the  PRFs  so  that  the  range  of  PRFs  required  is  less  extreme. 
We  will  illustrate  these  design  procedures  by  redesigning  the  baseline  sys¬ 
tems  under  the  assumption  that  four  isolated  parallel  channels  are  available 
at  frequencies  that  are  spaced  about  3%  apart.  (The  relationship  between  the 
new  adjusted  PRF  and  the  original  PRF  is,  for  the  i-th  burst. 


new 

ri 


X  f 
o  ri 


ong 


/X. 

1 


where  X  >  X . . ) 
o  1 

3.1.2  Design  for  Rotating  1-D  Phased  Array 

We  may  redesign  the  waveforms  for  the  long-range  zone  by  merely  multi¬ 
plying  the  parameter  6  by  q ,  the  number  of  channels.  (As  long  as  the  resulting 
number  of  bursts,  p,  is  an  integer  multiple  of  q  then  the  analysis  is  valid.) 
Thus  we  have 


721.6 

for  p=4. 

q=4 

1443.2 

for  p=8. 

Si 

II 

p  =  .005 

The  impact  on  the  design  is  phenomenal.  We  find  for  6=721.6  that  K=43,  p=-7 
is  optimum  and  for  6=1443.2  that  K=57,  p=8  is  optimum.  Since  we  are  constrained 
to  p=8  for  efficiency,  8  must  lie  in  the  range  800  to  1950,  and  K  must  be  45  to 
57.  If  we  assume  K  is  50  and  that  the  basic  PRF  is  750  Hz  at  the  longest  wave¬ 
length  then  we  find  that  this  wavelength  must  be 

X  =  2v  /f  K  =  5.3  cm 
max  ro 


which  is  much  less  than  half  the  admissible  wavelength  for  a  singl e -channel 
system.  The  frame  time  is,  however,  16  seconds  for  the  long-range  search! 
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It  is  more  reasonable  to  restrict  the  number  of  bursts  to  4  so  tli.it  the 
frame  time  is  8  seconds  or  less.  In  this  case  f  is  7111.6  and  tin  optimum  p 
is  7,  but  we  are  constrained  to  the  suboptimum  choice  of  p=4.  The  achievable 
improvement  factor  K  is  therefore  reduced  from  50  to  about  34.  The  corres¬ 
ponding  wavelength  i6 


X=2v  / f  K  =  7.8  cm 
max  ro 

which  is  at  least  40%  shorter  than  the  wavelength  that  results  for  the  one- 
channel  radar. 

For  the  near-range  mode,  we  have  the  parameters 

8  =  541.2 

p  =  0.02 

for  the  four-channel  system.  Again  restricting  the  number  of  bursts  to  p=4 
we  find  that  the  achievable  improvement  factor  is  14.  In  this  mode,  however, 
we  may  choose  p=8  without  an  increase  in  search  dwell  time,  and  the  more 
nearly  optimum  solution  of  K=17  results.  For  the  1875  Hz  basic  PRF  of  the 

near-range  mode,  the  achievable  wavelength  is  6.3  cm,  which  is  shorter  than 

the  wavelength  allowable  in  the  far-range  zone.  Thus  we  must  uvnel ruit.  our¬ 
selves  to  the  latter  wavelength  (X  =  7.8  cm)  and  work  backwards.  We  find  that 
the  result  ingK  factor  is  13.6.  As  with  the  near-range  designs  for  the  conven¬ 
tional  approaches  in  Section  2,  we  must  raise  K  to  the  next  integer,  which  is 
K=14.  The  basic  PRF  then  becomes  1832  Hz  and  the  near-range  search  is  ex¬ 
tended  to  82  km.  Since  K=14  can  be  achieved  with  p=4  this  solution  is  ac¬ 

ceptable.  It  leads  to  h  =  525,  which  is  97%  ol  the  value  used  when 

T  =1  sec  (i.e.,  T  =  0.97  sec), 
s  s 

3.1.3  Design  for  Stat ionary  2-D  Phased  Array 

Wc  may  redesign  the  system  with  the  fully-phased  array  in  an  analogous 
manner.  If  we  retain  the  allowance  of  1.5  sec  (out  of  8  sec)  for  the  long- 
range  search,  then  we  merely  need  to  mult  iply  tin.*  parameter  ■  by  the  number 
of  channels,  q=4.  The  design  parameters  are  therefore 


(•  1082.2 
p  =  0.005 


3-5 


resulting  in  the  optimum  solution  10=51,  p=7.  Since  p  must  be'  an  integer  mul¬ 
tiple  of  q  we  will  use  p=8  rather  than  p=7.  The  corresponding  K  is  still  51, 
however,  so  p=8,  K=51  is  the  solution.  The  resulting  wavelength  is  X  =  5.2  cm, 
which  is  35%  lower  than  the  mechanically  rotating  ease. 

If  we  now  use  this  wavelength  as  a  constraint  for  the  near-range  search, 
then  the  required  improvement  factor  is  K=20.5  for  a  1875  Hz  basic  PRF.  We 
may  therefore  relax  the  basic  PRF  to  1832  Hz  so  that  K  is  21  and  the  maximum 
range  is  now  82  km.  The  corresponding  8  for  p=8  and  K=21  is  1608.  This  solu¬ 
tion  implies  that  2.2  seconds  are  needed  for  short-range  search,  which  is  more 
than  the  allowed  2-second  frame  time.  This  solution  is  therefore  clearly  in- 
admissable! 

The  difficulty  is  that  the  near-range  search  now  controls  the  choice  of 
wavelength,  as  we  shall  see.  We  choose  Tg=.638  sec  for  the  short-range  search 
at  each  antenna.  This  is  the  time  used  for  the  short-range  search  and  2-D 
phased  array  when  only  one  channel  can  be  used.  The  corresponding  8  was  115. 
Since  we  have  four  channels  this  is  multiplied  by  four  so  that  in  the  near¬ 
range  design  we  have 

8  =  460 

p  =  0.02 

These  parameters  lead  to  K=16  at  p=6.  For  the  allowable  solutions  we  have 
K=16,  p=8  and  K=13,  p=4.  Using  p=8  we  obtain  X=  6.7  cm.  This  is  longer  than 
the  value  X  *  5.2  cm  found  for  the  long-range  search.  Thus  using  the  same 
search  times  as  we  did  for  the  one-channel  case  we  find  that  the  short-range 
wavelength  must  be  used. 

If  we  use  p=4  (which  was  used  in  the  mechanical,  four-channel  design)  we 
find  that  the  short-range  case  now  yields  X=  8.2  cm  while  the  long-range  case 
yields  X  =  7.2  cm.  Again  the  short-range  wavelength  of  8.2  cm  must  be  used. 

But  this  is  larger  than  the  wavelength  that  was  allowed  with  the  mechanically 
scanned  antenna  using  4  channels  and  p=4.  That  is,  we  end  up  with  a  worse 
solution  with  the  better  antenna! 

This  difficulty  can  be  overcome  by  using  less  time  in  the  long-range 
search  and  thereby  providing  more  time  in  the  short-range  search.  Thus  wo 
shall  use  1  sec  (out  ot  8  sec)  for  the  long-range  search.  For  the  long- range 
searr h  we  now  have 


If 
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B  =  721.8 
p  =  0.005 

resulting  in  the  optimum  solution  K=43,  p=7 .  For  the  allowable  values  of  p=4 
and  p=8,  we  have  K=34  and  K=43,  respectively. 

For  p=8  the  resultant  wavelength  is  6.2  cm  which  is  about  50%  of  the  value 
which  can  be  achieved  with  one  channel.  For  the  wavelength  of  6.2  cm  we  find 
in  the  short-range  search  that  we  must  use  a  K=18.  For  p=8  the  resultant  8 
is  625,  which  corresponds  to  a  scan  time  of  0.85  sec  (out  of  2  sec)  for  the 
short-range  case.  We  therefore  devote  42.5%  of  the  available  time  to  short- 
range  search  and  12.5%  of  the  time  to  long-range  search,  for  a  total  of  55% 
of  the  time  being  devoted  to  search  modes  altogether. 

In  sum  by  using  a  four-channel  system  we  have  reduced  the  wavelength  from 
11.6  cm  to  6.2  cm  by  devoting  an  equal  amount  of  dwell  time  to  both  the  long- 
range  and  short-range  searches.  The  tradeoff  is  not  all  advantageous,  how¬ 
ever,  because  we  now  use  about  twice  as  many  bursts  for  blind-speed  resolution. 
This  means  that  the  average  power  required  for  detection  is  doubled,  which 
could  exceed  practical  constraints. 

Suppose  we  instead  constrain  the  number  of  bursts  to  p=4  in  light  of  the 
possible  power  problem.  Now  we  can  design  the  four-channel  system  using  one 
burst  per  channel.  Under  this  constraint  the  long-range  mode  requires  K=34. 

The  allowable  wavelength  for  the  basic  PRF  of  750  Hz  is  7.8  cm.  Thus  we  still 
have  a  very  worthwhile  improvement  over  the  11.6  cm  of  the  one-channel  case. 

Using  7.8  cm  as  the  wavelength  constraint  on  the  near-range  search,  we 
obtain  K=14,  and  a  PRF  of  f  =  1832  Hz.  The  short-range  search  is  now  ex¬ 
tended  to  R  =82  km.  The  corresponding  8  for  p=0.02  and  p=4  is  8=525. 
max  r 

Using  the  fact  that  we  are  using  four  channels  this  corresponds  to  a  search 
time  of  0.72  seconds  for  the  short-range  case.  Now  we  are  using  36%  of  the 
time  for  short- -ange  search  and  12.5%  of  the  time  for  long-range  search.  This 
means  we  are  using  48.5%  of  the  time  for  search  which  leaves  51.5%  of  the  time 
for  tracking.  This  is  an  improvement  over  the  49.35%  of  the  time  left  for 
tracking  in  the  one-channel,  while  the  wavelength  of  7.8  cm  is  33%  smaller 
than  the  one-channel  case. 
hfdule  of  Biwg  ts 

The  long-range  search  uses  the  following  PRFs : 
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f  ,  =  750.00  Hz 
rl 

f  „  =  733.82  Hz 
r2 


r3 


692.75  Hz 


f  .  =  569.53  Hz 
rA 

A  single  blind  zone  of  19. A 
of  pulses  at  each  PRF  are 


Hz  width  remains  with  this  set  of  PRFs. 


The  number 


n^  =  11 

n ^  =  10 

n .  =  8 
A 

The  short-range  search  uses  the  following  PRFs 

f  ,  =1832  Hz 
r  1 

f  =  17AA.76  Hz 
r2 

f  ,  =  1582.75  Hz 
r3 

f  ,  =  1258.72  Hz 
rA 

A  single  blind  zone  of  37.  A  Hz  width  remains  with  this  set. 


3.1. A  Implementation 

We  have  seen  that  significant  performance  improvement  can  be  achieved  by 
using  a  multi-channel  radar  system  in  which  the  bursts  for  blind-speed  resolu¬ 
tion  are  implemented  in  parallel  on  different  RF  carriers.  We  have  not  addressed, 
however,  the  difficulty  of  implementing  such  a  system. 

Two  difficulties  are  evident.  First,  a  multi-channel  system  evidently  re¬ 
quires  a  duplication  of  RF  equipment  such  that  all  channels  can  be  used  at  the 
same  time.  This  requirement  may  be  difficult  to  meet  with  a  mobile  system. 
Moreover,  the  equipment  complexity  may  be  excessive.  Second,  there  is  a  severe 
problem  of  isolation  among  the  separate  channels.  Since  different  PRFs  are 
used  in  the  different  channels,  a  lack  of  adequate  isolation  could  easily  lead 
to  eclipsing  problems.  That  is,  reception  on  one  channel  may  be  required  at 
the  same  time  as  transmission  on  another  channel.  The  only  source  of  isolation 
is  the  separation  in  carrier  frequency  among  the  channels.  Our  experience  indi¬ 
cates  that  adequate  isolation  is  not  achieved  by  frequency  alone.  Instead,  it 
is  usually  necessary  to  use  separate  antennas  for  transmission  and  reception. 
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Since  weight  and  size  are  already  problems  for  this  mobile  radar,  such  a  solu¬ 
tion  is  probably  not  admissable. 

In  conclusion,  unless  the  technology  of  the  radar  equipment  can  be  ad¬ 
vanced  far  beyond  its  present  state,  particularly  in  the  area  of  transmit/ 
receive  isolation,  the  multi-channel  approach  outlined  here  will  be  impractical. 

3.2  SIMULTANEOUS  BLIND-RANGE/ BLIND-SPEED  RESOLUTION 
3 ■ 2 . 1  Background 

In  the  past  there  has  existed  a  consensus  that  ground-based  radars  which 

must  operate  in  a  heavy  clutter  environment  must  be  designed  without  range 

ambiguities  within  the  maximum  detection  range.  In  the  presence  of  range 

ambiguities,  a  target  to  be  detected  near  the  maximum  detection  range  would 

have  to  compete  against  clutter  folded  over  from  areas  very  close  to  the  radar. 

3  2 

The  K  law  for  ground  clutter,  and  in  the  case  of  rain  the  R  law  for  rain 
clutter,  would  enhance  the  clutter  relative  to  the  signal  to  such  a  degree 
that  clutter  suppression  would  be  impossible.  Thus  it  has  become  common 
practice  to  design  ground-based  radars  without  range  ambiguities. 

On  the  other  hand,  the  consequence  of  a  range-unambiguous  design  is  the 
apnearance  of  a  large  number  of  blind-speed  regions  within  which  aircraft  and 
missiles  cannot  be  detected.  Since  targets  may  go  at  constant  range  rate  for 
relatively  long  times,  in  particular  if  they  are  headed  toward  the  radar, 
blind  speeds  are  unacceptable  for  a  high-performance  radar.  Thus  one  must 
introduce  schemes  for  avoiding  blind  speeds,  which  means  multiple  PRFs.  If 
this  is  done  inefficiently,  one  is  wasting  such  important  radar  resources  as 
average  transmitter  power  and  dwell  time  on  the  target.  In  particular,  with 
respect  to  dwell  time,  it  is  found  that  avoiding  blind  speeds  may  increase 
the  required  dwell  time  to  sueh  an  extent  that  one  is  forced  to  choose  a  low 
cat  t  ier  frequency,  such  as  I, -band  or  even  liHF.  These  radars  then  have  poor 
angular  tracking  accuracy,  they  are  subject  to  severe  multipath  effects,  they 
are  more  readily  jammed,  and  they  have  an  insignificant  capability  for  target 
identification. 

It  appears  that  the  optimum  solution  (that  is,  the  solution  that  allows 

the  highest  carrier  frequency)  is  to  allow  both  blind  ranges  and  blind  speeds 

and  then  provide  sets  of  PRFs  that  allow  the  resolution  of  both  types  of 

3 

ambiguous  clutter.  The  problem  of  the  R  -enhanced  st length  of  range-ambiguous 
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clutter  presents  a  new  requirement  for  low  sidelobes  which  may  not  always  be 
achievable  in  practice.  Moreover,  the  elimination  of  all  blind-speed  and 
blind-range  zones  may  not  be  possible  in  general.  But,  as  we  shall  see,  the 
performance  potentials  are  so  high  from  a  combined  system  that  it  warrants 
our  serious  consideration  at  the  least  . 

3.2,2  The  Basic  Principle 

The  conventional  approach  to  PRF  selection  in  a  heavy  clutter  environ¬ 
ment  for  a  ground-based  radar  is  to  choose  the  fundamental  PRF  so  that  the 
first  range  ambiguity  will  be  at  the  maximum  detection  range.  If  multiple 
PRFs  are  required  to  avoid  blind  speeds,  then  all  other  PRFs  will  be  lower 
than  the  fundamental  one  so  that  there  will  be  no  range  ambiguities  closer 
than  the  maximum  detection  range.  The  wisdom  behind  this  choice  is  that 
range  ambiguous  operation  causes  an  enhancement  of  the  clutter,  which  is  to 
be  avoided. 

However,  in  spite  of  the  enhancement  of  clutter  there  is  much  to  be 
gained  by  range-ambiguous  operation.  The  additional  clutter  could  conceivably 
be  suppressed  by  Doppler  filtering,  and  the  higher  PRFs  make  the  blind-speed 
avoidance  problem  much  easier. 

In  order  to  illustrate  how  clutter  is  enhanced  by  range-ambiguous  oper¬ 
ation,  let  us  assume  that  the  fundamental  PRF  is  5  kHz,  so  that  the  first 

range  ambiguity  will  be  at  30  km.  In  volumetric  clutter  the  received  power 

2 

from  clutter  will  be  proportional  to  (1/R)  for  the  first  30  km.  At  R=35  km 

2 

the  received  power  will  be  proportional  to  (1/35)  plus  the  ambiguous  compo¬ 
nent  at  R=5  km,  so  that  the  total  clutter  power  will  be  proportional  to 
2  2 

(1/35)  +  (1/5)  .  In  general,  we  can  define 

R  -  target  range 

R^  =  range  of  first  ambiguity 

N  =  largest  integer  less  than  (R/R^) 

R'  =  R  -  NR. 

A 

The  received  clutter  power  will  then  be  proport ional  to 


If  we  compare  the  received  clutter  power  with  the  component  at  range  R,  we  can 
define  an  enhancement  factor 


which  is  illustrated  in  Figure  3-1. 

The  Doppler  filtering  processor  must  be  designed  to  suppress  the  enhanced 
clutter.  For  example,  if  the  processor  is  capable  of  handling  an  additional 
20  dB  of  clutter,  we  can  read  from  Figure  3-1  and  determine  that  the  particu¬ 
lar  pulse  burst  will  be  blind  at: 

R  <  R  <  1.10  R 
A  A 

2R,  <  R  <  2.25  Ra 

A  A 

3RA  <  R  <  3.50  Ra 

A  A 

4R  1  R  5  4.70  Ra 

A  A 

If  the  processor  could  handle  an  additional  30  dB  of  clutter  (10  more  than 
above)  we  can  limit  the  blind-range  interval  to  about  10%  of  R  up  to  the 

A 

fourth  range  ambiguity.  In  Figure  3-2  we  show  the  combined  range-Doppler 
blind  zones  for  a  2-PRF  burst  of  7  and  5  kHz  with  a  blind-range  interval  of 
0.10  R^  and  a  bl i nd-Doppler  interval  of  2  kHz.  Note  that  the  first  range  at 
which  both  PKFs  are  blind  is  7  times  the  first  blind-range  of  the  fundamen¬ 
tal  FRF  (7  kHz),  which  is  150  km.  The  first  Doppler  at  which  both  PRFs  are 
blind  is  14  kHz.  There  are  also  local  zones  at  which  both  PRFs  will  be 
blind,  but  these  zones  amount  to  only  about  5%  of  the  total  range-Doppler 
area  of  interest.  Additional  PRFs  can  be  used  to  improve  coverage,  or  to 
reduce  the  clutter  enhancement. 

We  have  expended  a  great  deal  of  effort  in  attempting  to  develop  syste¬ 
matic  procedures  based  on  these  principles,  but  the  choice  of  PRFs  for  oper¬ 
ation  in  range-ambiguous  clutter  is  not  yet  solved.  There  are  many  consid- 
erat  ions  such  as  the  maximum  detection  range,  first  blind  speed,  Doppler 
width  of  clutter,  and  clutter  suppression  capabilities  via  Doppler  filtering. 
We  lie  I  l eve  that  a  systematic  approach  could  be  evolved,  and  the  result  would 
be  superior  to  the  conventional  approach  of  unambiguous  operation.  We  will 
demonstrate  the  possible  improvements  through  the  following  intuitive  look  at 
the  approach. 


3-13 


3.2,3  Achievable  Performance  Improvement 

The  given  quantities  in  any  such  system  are  the  maximum  detection  range, 

R  ,  the  maximum  target  speed  to  be  accommodated,  v  ,  and  the  maximum 
max  max 

range  rate  spread  of  the  clutter,  AR^.  We  assume  that  the  ambiguous  range 
is  to  be  some  fraction  of  the  maximum  detection  range. 


R  =  R  /A 
a  max 


(3-1) 


where  it  would  be  desirable  for  A  to  be  no  larger  than  unity.  The  case  to 
be  discussed  is  where  A  must  of  necessity  be  larger  than  unity. 

From  (3-1),  the  repetition  frequency  is 

f  =  c/2R  =  cA/2R  ,  (3-2) 

r  a  max 


and  by  multiplication  by  A / 2  we  obtain  the  ambiguous  range  rate 

R  =  cAX/4R  .  (3-3) 

a  max 

The  normalized  basic  PRF  of  the  system  is  simply  the  ratio  of  R  and  the 

3 

range  rate  spread  of  the  clutter,  AR^,  so  that  from  (3-3)  we  have 


x 


R  /AR  „ 
a  cl 


cAX/4AR 


R 

ci  max 


(3-4 ) 


It  is  seen  from  (3-4)  that  the  normalized  PRF  decreases  with  the 
wavelength.  Since  blind-speed  resolution  is  impossible  when  x  is  much 
smaller  than  two,  a  decrease  of  the  wavelength  must  beyond  some  point  be 
offset  by  an  increase  in  the  value  of  A,  since  the  other  parameters  in  (3-4) 
are  fixed.  A  value  of  A  exceeding  unity,  however,  implies  the  introduction 
of  range  ambiguities,  so  that  the  system  of  multiple  PRFs  must  be  designed 
to  avoid  both  blind  speeds  and  blind  ranges.  The  interesting  question  is 
whether  there  are  any  limits  on  the  maximum  value  of  A,  and  what  these  limits 
are  if  they  do  exist.  In  other  words,  how  do  we  find  sets  of  PRFs  for  com¬ 
bined  blind-speed  and  blind-range  resolution  which  allow  us  to  achieve  high 
values  of  A,  and  are  there  any  fundamental  restrictions? 

On  the  basis  of  a  limited  study  of  the  problem,  we  have  come  to  the 
following  conclusions.  Although  there  is  a  systematic  way  of  designing  a 
set  of  PRFs  for  complete  elimination  of  blind  speeds,  this  is  not  possible 
for  the  two-dimensional  case  of  combined  blind-speed  and  blind-range 
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elimination.  Blind  regions  in  the  range/range-rate  plane  can  be  eliminated, 
or  made  insignificantly  small,  for  close  ranges.  With  increasing  range  an 
increasing  percentage  of  the  range/range-rate  plane  will  be  blind,  until 
beyond  some  range  essentially  the  entire  region  is  blind.  Thus  there  appears 
to  be  a  fundamental  restriction  on  bow  large  a  range  interval  can  be  accom¬ 
modated  with  acceptable  blind  regions. 

It  turns  out  that,  with  a  good  design  of  the  set  of  PRFs,  the  range  in¬ 
terval  over  which  blind  regions  can  be  reduced  to  acceptable  level  is  in¬ 
versely  proportional  to  the  range-rate  spread  of  the  clutter.  The  larger  the 
range-rate  spread,  the  smaller  the  range  interval  that  can  be  "cleared". 

This  maximum  operational  range  can  be  written  as 


R 

op 


cA 


4AR 


cSL 


(3-5) 


where  the  factor  n  depends  on  the  design  of  the  PRF  set. 

We  have  found  in  the  limited  research  that  a  factor  n=l  is  easy  to 
achieve,  but  that  the  blind  regions  increase  rather  rapidly  in  relative  size 
when  the  value  of  unity  is  exceeded.  There  appear  to  be  no  problems  of  find¬ 
ing  PRF  sets  for  n=l-2,  and  perhaps  n=1.3.  It  is  more  difficult,  but  possible, 
to  find  sets  with  n=1.5.  The  general  trend  is  that  larger  values  of  n  can  be 
achieved  with  a  larger  number  of  PRFs,  evidently  at  the  price  of  increased 
system  complexity  and  power  inefficiency.  We  shall  postpone  a  discussion  of 
how  the  PRF  sets  are  found,  and  simply  assume  here  that  the  value  of  n  that 
can  be  realized  is  limited.  Although  our  research  is  not  broad  enough  to 
permit  a  completely  reliable  estimate  for  the  maximum  value  of  n,  we  shall 
here  assume  as  a  reasonable  estimate 


1.5 


(3-6) 


We  return  to  our  system  design  problem,  which  means  finding  the  value 
of  the  normalized  PRF  x  of  (3-4).  Since  x  must  bo  larger  than  two  but  a 
small  wavelength  is  desired,  we  must  increase  the  value  of  A.  Thus  we  con¬ 
sider  next  the  question  of  how  large  A  can  be  made. 

We  evidently  must  make  the  operational  range  K  of  (3-5)  equal  to  the 


maximum  detection  range  R 


Substituting  R  for  K  in  (3-4)  gives 
op  max 
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so  that 


x  =  A/n 


(3-7) 


A  =  x  n 
max  max 


(3-8) 


If  we  substitute  A  into  (3-4),  we  find  t lie  luimiraum  wavelength  lor  whirl) 
max  b 

combined  blind-range  and  blind-range  rate  resolution  is  possible. 


A  .  =  4AR  R  /cr\ 

mm  ct  max  max 


(3-9) 


With  n  =  1.5,  we  have 
max 

A  .  =  2.7  AR  „  R  /c 

mm  cl  max 


(3-10) 


As  an  example,  suppose  AR  „  =  20  m/sec  and  R  =  200  km.  The  minimum 

cl  max 

wavelength  for  which  the  problem  of  blind  ranges  and  blind  speeds  can  be 

overcome  then  is  A  .  =3.6  cm.  Since  the  minimum  value  of  x  is  two,  we  see 

mm 

from  (3-8)  that  such  a  design  would  have  at  least  a  triple  range  ambiguity. 
In  practice,  the  wavelength  may  be  fixed.  If  it  is  not  fixed,  then  we 


may  estimate  from  (3-10)  the  minimum  usable  wavelength  from  the  point  of  view 
of  avoiding  clutter  problems,  and  assess  other  problems  at  this  wavelength 
(such  as  rain  clutter  and  rain  attenuation  strengths).  This  will  then  lead 
to  a  decision  about  the  wavelength  to  be  used. 

Once  the  wavelength  is  chosen,  we  can  from  (3-4)  find  the  relation  be¬ 
tween  x  and  A.  A  must  be  large  enough  to  yield  a  value  of  x  exceeding  two, 
and  the  larger  the  value  of  x,  the  fewer  the  number  of  PRFs  that  must  be 
used.  Thus  we  want  A  Large.  On  the  other  hand,  the  more  closely  A  approaches 
its  maximum  value  as  given  by  (3-8),  the  more  serious  the  residual  clutter 

problems  will  be.  Thus  the  value  of  A  is  bracketed  between  A  .  eorrespond- 

nun 

ing  to  a  value  of  x=2,  which  from  (3-4)  is 


A  .  =  8aR  R  /cA 

min  cE.  max 


(3-11) 


and  the  value  A  of  (3-8).  However,  the  latter  constraint  means  only  that 
max 

x  must  be  made  large  enough  to  accommodate  the  chosen  value  of  A.  A  high 
value  of  A  means  that  there  are  many  effective  range  ambiguities.  The  sim¬ 
plified  derivation  on  which  the  above  conclusions  are  based  does  not  fully 
take  into  account  the  increase  In  clutter  problems  with  increasing  order  of 
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the  range  ambiguity.  Thus  we  want  to  keep  A,  which  means  the  number  of 

range  ambiguities,  as  small  as  possible;  that  is,  close  to  the  value  of 

A  .  of  (3-11). 
min 

In  the  design  procedure,  we  find  from  (3-11)  the  value  of  A  .  ,  and 

min 

from  (3-7)  the  corresponding  value  of  x.  We  decide  at  this  point  how  criti¬ 
cal  the  clutter  problem  and  the  situation  are,  by  deciding  whether  we  want 
to  design  for  n=1.2,  1.2,  or  a  higher  value.  The  more  PRKs  that  are  used,  the 
higher  the  admissible  value  of  n,  but  also  the  higher  the  cost.  We  survey 
the  available  solutions  for  the  sets  of  PRKs,  and  decide  on  the  best  compro¬ 
mise  between  x  and  p.  What  is  still  needed  is  the  set  of  "best"  solutions 
for  x  and  the  associated  other  normalized  PRF  values. 

It  appears  that  the  design  of  an  optimum  set  of  PRKs  is  not  possible. 

We  have  found  one  set  of  solutions  by  a  combination  of  reasoning  and  computer 
trials.  We  have  found  other  solutions  by  starting  with  the  sets  of  PRKs  used 
for  resolution  of  blind  speeds  only,  and  multiplying  the  normalized  PRKs  by 
some  factor.  This  has  revealed  the  general  trend  of  results  as  discussed 
earlier,  yet  some  solutions  are  better  than  others  and  locally  reverse  the 
trend.  The  best  procedure  for  a  generalized  system  of  solutions  appears  to 
be  t Ho  following:  Start  with  the  blind-speed  only  sets  of  PRKs  for  three, 
four,  and  five  PRKs,  and  in  each  case  for  different  values  of  the  extension 
factor  for  the  blind  speed.  Also  multiply  these  normalized  PRKs  by  some 
selected  factors.  For  these  diverse  sets  of  PRFs  we  then  find  the  value  of 
r,  lor  which  percentage  blindness  at  the  largest  range  exceeds  a  certain 
threshold.  We  then  select  those  solutions  for  which  ti  is  largest  for  a  given 
number  of  PRKs  and  a  given  extension  factor  (or  equivalent  measure).  The  re¬ 
sulting  table  of  solutions  is  the  one  to  be  used  in  the  system  design.  The 
.same  solutions  are  used  in  an  adaptive  system,  where  the  clutter  conditions 
are  changing  and  range  ambiguous  operation  is  necessary  under  some  conditions 
and  in  some  beam  positions. 

3.2.4  Application  to  Radar  Design  Concepts 

The  difficulty  of  finding  a  systematic  solution  to  the  combined  blind- 
speed/blind  -  range  problem  has  led  ns  to  t  lie  conclusion  that  a  further  study 
is  beyond  the  scope  of  the  present  program.  However,  We  have  recently  had 
the  opportunity  to  study  the  problem  of  a  1 ow-probabi 1  it y-of~ intercept  ( HP  I ) 


t 


3-17 


radar  In  which  the  principles  just  discussed  are  applicable.  Our  study  of 
the  LP1  problem  has  led  to  the  conclusion  that  the  combined  blind-range/ 
blind-speed  problem  must  be  solved  if  future  radars  are  to  be  made  signifi¬ 
cantly  less  vulnerable  to  anti-radiation  missiles.  We  therefore  review  the 
tradeoffs  here  to  demonstrate  that  a  separate  program  is  needed. 

Background 

Our  study  of  the  LPI  problem  has  shown  that  a  common  assumption  made  in 
these  types  of  investigations  is  that  the  ARM  will  use  essentially  the  same 
type  of  receiver  against  LPI  radars  as  is  being  used  against  existing  radars. 
Specifically,  the  assumption  is  made  that  the  ARM  will  not  use  a  high  degree 
of  noncoherent  integration  of  the  radar  signal.  We  believe  that  this  is  an 
entirely  unrealistic  assumption.  Since  the  signal  of  a  conventional  radar 
can  be  detected  with  extreme  ease  by  an  intercept  receiver,  there  is  no  justi¬ 
fication  for  designing  any  but  the  simplest  type  of  intercept  receiver. 
However,  with  LPI  developments  under  way,  it  is  reasonable  to  assume  that 
better  intercept  receivers  will  be  designed.  In  particular,  it  is  simple  and 
inexpensive  to  design  the  intercept  receiver  so  that  it  uses  long-term  non¬ 
coherent  integration.  However,  with  current  radar  designs  and  conventional 
design  approaches  one  can  show  it  to  be  impossible  to  prevent  easy  homing  on 
part  of  the  ARM.  Although  it  is  true  that  one  can  deploy  radar  decoys  in 
order  to  protect  the  radar,  these  decoys  become  rather  expensive  and  vulner¬ 
able  if  designed  for  conventional  radars. 

At  MARK  Resources  we  have  studied  the  problem  and  have  found  that  LPI 
performance  can  be  obtained  if  nonconventional  design  approaches  are  being 
used  for  the  radars.  In  fact,  these  approaches  appear  to  be  so  effective 
that  it  may  be  unjustified  to  design  modern  radars  in  accordance  with  the  old 
methods,  including  the  ones  we  have  developed  at  MARK  Resources.  Hence  it 
appears  to  us  extremely  important  to  extend  the  present  problem  on  radar 
adaptivity  to  the  inclusion  of  LPI.  We  see  a  way  of  modifying  our  design 
methods  to  where  effective  LPI  radars  can  be  designed  just  as  efficiently  and 
systematically  as  we  have  done  previously  for  "ordinary"  radars. 

Brief  Description  of  the  New  Approach  to  Radar  Design 

There  are  several  departures  from  conventional  radar  design  which,  taken 
together,  allow  the  implementation  of  a  true  LPI  radar.  By  this  we  mean  a 
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radar  that  prevents  successful  homing  by  an  ARM.  We  shul  1  brieiiy  des.  ri!u 
these  new  design  features. 

As  perhaps  the  most  important  features,  the  radar  must  be  designed  so 
that  it  operates  intermittently.  For  example,  it  might  search  for  half  a 
second  and  not  transmit  at  all  for  the  next  half  second.  It  is,  of  course, 
obvious  that  such  intermittent  operation  is  highly  desirable.  What  is  not 
so  obvious  is  how  to  perform  all  the  required  radar  functions  within  the  re¬ 
duced  time,  where  even  with  conventional  designs  it  is  difficult  to  provide 
sufficient  time  for  all  functions. 

Our  study  of  the  problem  has  shown  that  dwell  time  requirements  can  be 
drastically  reduced  by  going  to  range-ambiguous  operation.  The  ambiguous 
range  is  chosen  much  smaller  than  the  maximum  detection  range.  We  have  found 
that  the  unsolvability  of  the  clutter  enhancement  problem  due  to  the  ambigu¬ 
ous  foldover  is  a  myth.  Although  penalties  must  be  paid,  these  penalties  can 
be  made  to  be  acceptable.  The  overall  result  is  that  the  dwell  time  require¬ 
ments  are  reduced  by  a  factor  roughly  in  the  order  of  the  multiplicity  of  the 
range  ambiguity.  In  other  words,  if  we  introduce  five  range  ambiguities 
within  the  maximum  detection  range,  we  would  expect  a  dwell  time  reduction 
by  a  factor  of  roughly  five. 

The  scheme  for  resolving  blind  speeds  for  range-unambiguous  operation 
must,  be  replaced  by  one  for  resolving  the  combination  ol  blind  speeds  and 
blind  ranges.  We  have  done  enough  work  to  know  that  the  problem  is  solvable, 
and  we  have  investigated  particular  approaches.  We  also  understand  the  limi¬ 
tations,  but  these  limitations  are  beyond  the  typical  practical  combinations 
of  maximum  target  speeds  and  maximum  detection  ranges.  What  is  not  available 
is  a  unified  L reory  of  optimum  algorithms  for  selecting  multiple  PRFs  for 
combined  blind-speed/blind-range  resolution  for  minimizing  dwell  time 
requi rements . 

A  second  important  feature  is  Lhe  use  ol  many  range  zones  during  the 
search.  The  farthest  range  zone,  for  example,  requires  the  highest  trans¬ 
mitter  power,  but  the  search  f rame  time  for  this  far  zone  can  be  very  long, 
perhaps  close  to  a  minute  for  really  Jong  search  ranges.  A  high  power 
signal  that  is  observed  only  once  every  minute  evidently  is  of  no  use  for 
an  ARM.  As  the  search  zones  get  closer  to  the  radar,  the  search  frame  t  ime 
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must  be  reduced.  However,  the  transmitter  power  will  be  reduced  in  accordance 
A 

with  the  R  law.  Thus,  as  the  frame  times  become  shorter,  the  signal  will 
be  more  and  more  difficult  to  detect.  We  may  use  a  minimum-range  search  zone 
for  which  the  largest  detection  range  is,  say,  only  20  km,  with  a  search 
frame  time  of  only  one  second.  The  radar  power  then  will  be  many  orders  of 
magnitude  lower  than  for  the  maximum  range;  so  low  in  fact,  that  successful 
homing  on  the  part  of  the  ARM  can  be  pi  evented. 

We  have  used  different  search  zones  in  our  earlier  designs,  but  for  dif¬ 
ferent  purposes  and  subject  to  different  constraints.  ror  LPI  radar,  for 
example,  we  want  to  time  the  search  over  all  zones  so  that  when  the  search 
for  different  zones  coincides,  it  can  be  done  with  a  single  transmission, 
just  as  with  an  ordinary  radar.  This  implies,  however,  that  the  multiple  PRF 
scheme  for  blind-speed/blind-range  resolution  is  appropriately  designed.  A 
consecutive  search  over  different  search  zones  must  be  avoided  in  order  to 
achieve  a  longer  dead  time  for  the  radar. 

Another  important  design  constraint  concerns  the  presence  of  multiple 
radars.  It  is  usually  argued,  with  much  hand-waving,  that  a  radar  with  only 
a  minimal  LPI  capability  might  prevent  homing  by  an  /ARM  since  there  will  be 
more  than  one  radar  in  the  battlefield.  A  closer  ex  imination  of  this  point 
shows  that  the  presence  of  multiple  radars  can  indeed  be  utilized  for  im¬ 
proved  LPI  performance,  but  only  if  this  is  taken  into  account  in  the  design 
procedures.  For  example,  use  of  a  very  high  carrier  frequency  will  allow 
the  ARM  to  implement  a  sharp  beam,  for  better  angular  resolution  oi  multiple 
radars.  On  the  other  hand,  much  of  this  advantage  can  be  offset  by  the  tre¬ 
mendous  differences  in  the  power  levels  used  by  these  radars  for  the  differ¬ 
ent  search  zones.  A  radar  in  the  outskirts  of  the  ARM  beam  and  transmitting 
at  the  power  level  for  the  far  zone  can  appear  much  stronger  for  the  .ARM 
than  a  radar  within  the  beam  and  transmitting  at  the  power  level  for  the 
close  zone.  The  resulting  phase  center  wander  is  another  contribution  to 
the  defeat  of  the  ARM, 

Recommendations 

These  preliminary  results  indicate  that  a  program  should  be  initiated 
for  extending  our  design  procedures  to  radar  witli  LPI  capability.  The  central 
issue  is  that  of  combined  blind-speed  and  blind-range  avoidance  by  multiple 


PRFs,  and  its  interrelation  with  clutter  enhancement  due  to  range  ambiguities 
Optimum  methods  of  dividing  the  range  interval  into  multiple  search  zones, 
with  a  particular  view  toward  the  compromise  between  reduction  of  the  dwell 
time  and  increased  clutter  enhancement,  must  be  developed.  The  interrelation 
with  search  frame  times  as  a  function  of  distance  from  the  radar  and  power 
levels  as  seen  by  the  ARM  must  be  studied,  and  optimum  design  procedures  must 
be  developed.  In  essence,  we  should  arrive  at  a  point  where  the  design  of  an 
LP1  radar  can  be  carried  out  in  the  same  systematic  fashion  as  for  a  conven¬ 
tional  system. 


4.0  DESIGNS  FOR  ENVIRONMENT -ADAPTIVE  RADAR 


The  designs  reviewed  thus  far  have  been  cognizant  of  the  environment 
but  have  not  been  adaptive  to  it.  For  example,  the  systems  have  been  designed 
to  handle  ground  clutter  to  some  maximum  range  (say,  80  or  100  km),  but  have 
not  been  designed  to  readjust  themselves  as  this  maximum  range  changes  with 
time  (or  with  time-varying  azimuth  angle).  The  final  big  step  in  performance 
improvement  is  expected  from  this  and  other  adaptive  readjustments  to  the  en¬ 
vironment.  Just  how  this  is  done  is  detailed  in  this  section  along  with  a 
preliminary  estimate  of  the  expected  performance  improvement. 

4.1  REAL-TIKE  ADAPTATION  OF  WAVEFORM  DESIGN 

4.1.1  System  Approach 

The  Need  for  Adaptivity 

We  have  previously  determined  that  the  solution  of  the  blind-speed  prob¬ 
lem  dominates  the  radar  design.  It  is  therefore  natural  that  a  readjustment 
of  the  solution  to  a  changing  environment  would  dominate  the  design  of  an  adap¬ 
tive  radar.  Since  the  blind-speed  problem  is  solved  through  proper  waveform 
design,  the  basic  problem  in  developing  an  adaptive  radar  therefore  lies  in 
developing  a  means  of  redetermining  and  readjusting  the  waveforms  as  the  en¬ 
vironment  changes. 

There  are  basically  three  environmental  characteristics  that  greatly  in¬ 
fluence  the  waveform  design:  range  extent  of  the  ground  clutter,  Doppler  spread 
of  the  clutter,  and  intensity  of  the  clutter  backscatter.  These  characteris¬ 
tics  are  clearly  time  varying. 

On  a  short-term  basis  they  are  time  varying  because  of  the  rapidly 
changing  azimuth  look  angle  of  the  radar  antenna.  The  search  is  performed 
with  a  beam  that  is  narrow  in  azimuth  and  scans  sequentially  through  the 
search  volume  (e.g. ,  through  360  degrees  of  azimuth).  Each  azimuth  angle 
will  generally  have  different  clutter  characteristics  so  that  I  he  waveform 
must  be  capable  of  rapid  readjustment  within  each  search  frame.  For  example, 
a  northerly  wind  will  cause  rain  to  present  a  high  mean  Doppler  and  Doppler 
spread  to  the  radar  when  the  radar  is  pointed  to  north  or  south,  but  a  low 
mean  Doppler  and  Doppler  spread  when  pointed  to  east  or  west.  As  a  eon.se- 
quence  the  total  clutter  spread,  including  the  ground  clutter  which  has  zero 
Doppler,  will  be  maximum  at  north  and  south  azimuths  and  minimum  at  east  and 
west  azimuths. 


The  clutter  characteristics  are  also  time-varying  on  a  long-term  basis. 
That  is,  over  several  search  frames  the  weather  may  cliange  enough  to  cause 
the  clutter  characteristics  to  change  even  at  a  f ixed  azimuth  angle.  High 
winds  may  increase  clutter  spread,  and  high  moisture  may  increase  the  clutter 
intensity  of  vegetation.  Moreover,  an  adversary  may  dispense  chaff,  which  will 
alter  the  required  waveforms  tremendously . 

Inplcrrumtaiion  of  Adaptivity 

Evidently  the  earth-referenced  clutter  conditions  clrange  relatively 
slowly  with  time  even  though  the  radar-referenced  conditions  change  rapidly 
because  of  a  rapidly  changing  azimuth  angle.  This  fact  suggests  that  the 
waveform  changes  with  azimuth  angle  should  be  implemented  with  a  rapid  ]ouk-u? 
table,  where  the  contents  of  the  table  are  modified  relatively  slowly  as  the 
earth-referenced  environment  changes.  In  other  words,  real-time  waveform  re¬ 
design  is  a  relatively  slow  procedure  because  of  a  slowly  changing  environment, 
whereas  real-time  waveform  resetting  is  a  rapid  process  that  occurs  many  times 
within  a  search  frame  but  is  based  on  waveforms  that  have  been  previously  de¬ 
signed  . 

This  two-pronged  approach  to  the  adaptive  radar  is  illustrated  in 

Figure  4-1.  An  auxiliary  sensor,  which  may  actually  be  integrated  with  the 

radar,  measures  the  environment  and  thereby  determines  the  design  parameters 

for  the  waveform  designs.  For  an  azimuth  angle  0  the  relevant  parameters  will 

be  p(0),  R  ,  (0),  and  k  (0)  ,  where  p  is  the  ratio  of  clutter  spread  tc  maximum 
amb  p 

target  speed  of  interest,  R  is  the  range  extent  of  the  ground  clutter,  and  k 

is  the  filter  broadening  required  to  provide  adequately  low  sidelobes  for  sup- 

pit  using  the  sensed  clutter.  A  360°  search  sector  may,  for  example,  be  divided 

into  16  sectors  of  22  l/2‘  each,  with  separately  measured  values  for  p,  R  , 

amb 

and  k  in  each  sector.  Whenever  the  measured  values  in  a  sector  change  sig¬ 
nificantly  from  those  that  resulted  in  the  current  waveform  set,  a  new  set  is 
computed  by  means  of  a  real-time  algoritlim  for  waveform  design.  This  new  set 
then  replaces  the  old  set  in  the  memory  for  that  sector  so  that  the  new  set 
will  be  used  the  next  time  the  coi responding  azimuth  is  passed. 

The  payoff  for  these  real-time  readjustments  is  two-fold.  First,  the 
dwel 1  lime  required  per  search  frame  is  minimized  for  each  environment  so 
that  t hi>  remaining  time  available  for  tracking  is  maximized.  Second,  bee. 
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Figure  4-1.  Architecture  for  Adaptive  Radar 
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the  time  is  used  so  efficiently  the  allowable  wavelength  for  a  given  perfor¬ 
mance  level  is  reduced  below  the  wavelength  allowable  in  a  non- adaptive  system. 
In  theory,  the  judicious  use  of  adaptivity  allows  one  to  design  for  the  average 
conditions  rather  than  the  worst-case  conditions  with  the  knowledge  that  the 
adaptive  networks  will  adjust  to  adverse  conditions  when  needed.  In  practice, 
the  radar  will  be  designed  such  that  the  search  frame  time  is  not  exceeded  by 
the  search  dwell  time  even  in  the  worst  case.  Then  as  conditions  become  more 
benign  the  waveforms  will  be  redesigned,  thereby  freeing  more  and  more  time 
for  tracking. 

In  the  following  we  will  develop  the  formalism  for  the  waveform  redesign. 
First  we  will  outline  how  the  optimum  design  parameters  are  determined  (e.g., 
search  time,  number  of  bursts,  basic  PRF),  then  we  will  demonstrate  a  real¬ 
time  algoritlim  for  deriving  the  full  PRF  schedule  from  these  parameters. 

4.1.2  Optimum  Design  Parameters 

The  adaptive  radar  is  controlled  by  spatio-temporal  changes  in  the  environ¬ 
ment  as  summarized  by  just  three  environmental  parameters: 

R  =  range  extent  of  ground  clutter 

#  cl 

AR^  =  range-rate  spread  of  total  clutter 
S^£  =  Doppler  sidelobe  level  required  for  clutter  suppression 

We  shall  differentiate  between  long-term  time  variations  and  short-term  time 
variations  as  described  previously  by  redesigning  the  waveforms  as  a  function  of 
Ra  (0,t) 

arc£  ie,t) 

Sc£  (0,t) 


where  0  is  ihe  azimuth  angle  and  t  is  time.  Although  azimuth  angle  is  a  strictly 
rapidly  varying  function  of  time  for  a  sequential  scan,  we  shall  treat  each 
azimuth  independently  and  thereby  treat  azimuth  as  a  constant.  The  time-varying 
azimuth  angle  will  then  be  handled  by  the  commutator  illustrated  in  Figure  4-1. 

Tn  other  words,  we  will  divide  the  search  volume  into  N  azimuth  sectors  and 
perform  Independent  waveform  design  for  each  sector  on  the  basis  of  the  param¬ 
eters 


ARct(Vo 

ScdV‘> 


fur  n=l  to  N. 
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In  developing  the  real-time  procedures  for  waveform  reoptimization  we 
must  first  translate  the  three  environmental  parameters  into  parameters  that 
are  recognizable  in  our  automatic  blind-speed  algorithms.  Recall  that  we 
have  two  search  zones:  long-range  zone  and  short-range  zone.  The  parameters 
that  are  measured  by  the  auxiliary  sensor  must  also  be  measured  separately  in 
each  zone.  The  parameters  to  be  measured  for  the  short-range  zone  are 


R  (0  ,t) 

.  a  n 


(0n,t) 

(6n,t) 


S 

S 


where  the  subscript  S  indicates  thg  short-range  zone.  The  parameters  for  the 
long-range  zone  are 


Sci(°n,C)L 


(The  maximum  operating  range,  R  ,  is  also  relevant,  but  it  is  fixed,  unlike 

the  maximum  range  for  the  short-range  zone,  R^.) 

The  input  parameters  for  the  automatic  design  algorithms  are  KCS^.t)  and 

p(0n>t),  which  must  be  defined  separately  for  the  long-  and  short-range  zones. 

The  algorithms  will  then  result  in  the  parameters  a(Bf,t)  and  p(0^,t).  The 

search  dwell  time  for  the  sector  is  then  determined  from  a(Q  ,t)  from  a 

n 

knowledge  of  k  (O^.t).  The  time  remaining  for  tracking  is  then  evident.  The 
required  translations  are  therefore,  for  a  given  system  wavelength. 
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and  K(0  ,t)T  is  fixed  by  R 

n  L  max 

The  translation  of  maximum  ranges  into  the  improvement  factor  K  is  based 
the  implicit  assumption  of  avoiding  range  ambiguities.  Thus  the  basic  PRF  is 


on 


c/2R  (8  ,t)  for  short-range  mode 
a  n 


c/2R  for  long-range  mode 
max  ° 


The  improvement  factor  is  therefore 


— 7 —  R  (0  »t)  for  short-range  mode 
cA  a  n 


— 7 —  R  for  long-range  mode 

cA  max 


The  translation  of  clutter  spread  into  the  parameter  p  is  direct: 


p(0  ,t)  =  AR  .(0  ,t)/v 
n  cZ  n  max 


with  an  S  subscript  for  the  short-range  zone  and  a  L  subscript  for  the  long- 
range  zone. 

The  translation  of  required  sidelobe  level  into  filter  broadening  fac¬ 

tor  is  much  less  direct.  The  precise  translation  depends  on  the  precise 
filter  synthesis  used.  When  deep  sidelobe  suppression  is  required  over  only 
a  narrow  fraction  of  the  sidelobe  region,  the  broadening  is  relatively  be¬ 
nign  and  is  near  unity.  However,  when  the  suppression  is  required  uni¬ 
formly  over  the  entire  sidelobe  region,  the  broadening  can  be  wide.  Moreover 
a  restriction  to  relatively  simple  filter  types  leads  to  much  worse  than  min¬ 
imum  broadening.  For  the  design  purposes  we  will  conservatively  assume  a 
simple  cosine-on-a-pedestal  class  of  filters.  The  broadening  can  then  be 
written  as 

k  (0  ,t)  =  0.886(1  +  0.636[-  .3409  +  /.04545 (-SJ1)  -  .47466]2} 

P  n  dB 

where  is  the  required  (negative)  sidelobe  level  in  dB,  or 


SdB  =  10  log10  ScT(en’° 

(again  with  an  S  subscript  for  the  short-range  zone  and  a  L  subscript  for 
the  long-range  zone). 
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4.1.3  The  Use  of  Subroutines  for  Determining  K,j>  in  Real  Time 

For  each  new  R  .  AR  „  and  S  .  that  is  sensed,  we  proceed  to  input  o  t  lie 
a  c£  c  C 

new  K  and  p  as  follows. 

4v  R 
ma  x  a 

K  “  cA  "" 

From  the  initial  design  v  and  A  are  fixed  so  tliat  we  determine  K  i  10m  the 

y 

new  R  alone.  In  general  the  resulting  K  will  not  be  an  integer.  Since  K 

a  cj 

must  be  an  integer  we  use  the  next  highest  integer  as  already  discussed.  Then 
we  have 


K  =  smallest  integer  no  smaller  than 


4v  K 

ma  x  a 
cA 


The  value  of  p  is  then  obtained  from  subroutine  OPB. 
we  have  inputs  K  and 


In  this  subroul  i n» 


p  =  AR  „/v 

c£  max 


Thus  p  is  determined  from  our  new  value  of  AR  „  .  The  subroutine  outputs  the 

cE 

optimum  value  of  p.  Thus  we  have  the  required  values  of  K  and  p. 

The  new  value  of  S  .  is  not  needed  in  order  to  determine  the  now  values 
cr 

of  K  and  p.  It  does  however,  determine  search  dwell  time  in  the  following  man¬ 
ner.  The  subroutine  OPB  outputs  the  minimum  B  along  with  the  corresponding  op¬ 
timum  value  of  p.  We  can  relate  k  to  S  ,  from  the  relations  in  Section  4.1.7. 

p  c£ 

The  dwell  time  required  for  search  is  then  determined  from 


k  w  9 

-P-®  -8  -  .B 

k  v  ' 
a  max 


4._1.4  New  Subroutines  for  Determining  PRF  Schedules  in  Real  Time 
Need  for  a  New  Approach 

Once  the  improvement  factor  K  and  number  of  bursts  p  have  been  deteimined 

we  still  have  the  problejn  of  determining  the  precise  set  of  PRFs.  In  the  past 

k 

we  have  relied  on  Johnson's  tables  for  defining  tin's  set.  However,  there  are  a 

* 

Johnson,  op  cit. 
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number  of  reasons  why  this  is  not  the  best  approach  for  incorporation  in  our 
environment-adaptive  radar. 

First,  a  table  of  all  sets  of  PRFs  from  Johnson's  tables  would  be  very 
unwieldy.  Numerous  sets  of  PRFs  would  have  to  be  stored  by  the  radar  when 
all  useful  pairs  K,  p  are  considered. 

Second,  the  tables  are  not  complete  in  that  not  all  possible  K,  p  pairs 

are  included.  The  reason  for  this  shortcoming  is  that  some  pairs  have  no 

solution  for  which  all  speeds  are  unblinded.  Johnson  has  tabulated  only 

those  sets  of  PRFs  for  which  there  are  absolutely  no  blind  speeds  from  0  to 

v  .  In  practice,  however,  we  would  allow  a  few  narrow  blind-speed  zones 
max 

because  such  narrow  zones  really  involve  only  the  skirts  of  the  clutter 
spectrum.  As  a  result  their  inclusion  would  degrade  detection  performance 
only  slightly  overall,  especially  considering  that  detection  probabilities 
no  greater  than  0.9  are  expected  anyway.  Had  Johnson  included  solutions 
with  a  few  narrow  unblinded  regions,  his  tables  would  have  been  complete  and 
useful  to  us.  On  the  other  hand,  since  he  maintained  a  stricter  definition 
of  the  blind  zone, we  find  that  we  must  develop  more  practical  solutions  of 
our  own. 

Third,  not  only  does  a  looser  definition  allow  completion  of  the  tables 
but  it  also  leads  to  more  practical  solutions  in  the  cases  where  Johnson  does 
have  solutions.  For  our  application  a  few  narrow  blind-speed  regions  would 
be  admissible  if  these  regions  occurred  only  at  low  speeds.  With  this  allow¬ 
ance  for  some  blind  zones  we  find  that  we  can  evolve  sets  of  PRFs  which  use 
fewer  bursts  and  a  narrower  range  of  PRFs,  thereby  increasing  the  practi¬ 
cality  of  the  solutions. 

H.vi  unp  lc  ,?s.  the  New  Approach 

before  mathematically  detailing  our  new  approach  to  finding  the  sots  of 
l’KFs,  let  us  outline  the  philosophy  and  give  an  example.  The  automatic  al¬ 
gorithm  for  computing  the  PRFs  picks  PRFs  in  decreasing  values  such  as  to 
unblind  the  higher  ambiguities  first.  Thus  if  we  do  not  use  the  lowest  PRFs, 
then  the  lower  ambiguities  may  be  partially  or  wholly  blinded.  It  is  of 
interest  in  some  situations  to  see  what  regions  are  blinded  when  the  lower 
PRFs  are  not  used.  Because  these  regions  correspond  to  low  velocity  rates, 
it  is  unlikely  that  threatening  targets  will  lie  in  these  regions.  If  this 
is  the  case  t hen  we  clearly  gain  by  using  fewer  PRFs. 
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As  an  example  we  consider  the  baseline  design  loi  the  mechanically  ro¬ 
tating  set  of  antennas  in  the  lar-range  zone.  By  using  our  new  algorithm 
(to  be  described  next)  lor  determining  the  I’KKs,  we  find  that  an  acceptable 
set  is 

f  =  750  Hz 
r  1 

f  ,  =  722.7b 
r  2 

l  ..  =  668.26 
r  i 

f  .  =  622.94 
r4 

This  solution  differs  from  the  set  that  would  be  iound  i rum  Johnson's  tables 
in  that  a  narrow  blind-speed  zone  remains.  This  zone  is,  however,  only 
5.25  Hz  wide,  which  is  certainly  small  by  comparison  with  the  15  kHz  Doppler 
spectrum  being  searched.  Moreover,  the  zone  is  centered  at  a  low  speed 
(about  50  m/s).  The  existence  of  this  residual  blind  zone  is  therefore  ol 
little  practical  consequence. 

Our  new  solution  also  has  the  property  that  the  omission  of  the  lowest 
PRF  in  the  schedule  introduces  blind  zones  at  the  lower  speeds  tirit.  This 
property  is  illustrated  in  Table  4-1.  Evidently  the  omission  of  each  suc¬ 
cessive  PRF,  starting  at  the  lowest  PRF,  merely  introduces  more  blind  ..ones 
at  the  lower  speeds.  A  good  tradeoff  therefore  exists  between  the  number  of 
bursts  used  and  the  size  of  residual  blinded  regions  at  low  speeds.  Thus, 
when  transmitter  power  or  dwell  time  is  at  a  premium  one  may  opt  to  omit  the 
PRF  of  lowest  frequency. 

Description  of  the  IRF  Selection  Algovit  hw 

We  can  now  describe  our  new  procedure  for  selecting  the  PRFs  once  the 
parameters  K  and  p  have  been  selected.  We  shall  choose  the  frequencies  so 
as  to  first  eliminate  the  higher  blind  speeds  as  these  are  the  most  important 
regions  in  practice.  Thus  if  after  using  the  p  frequencies  there  are  still 
blind  regions,  they  will  he  at  the  lowest  blind  speeds.  With  this  precedin'* 
we  can  also  examine  how  many  extra  blind  regions  are  created  at  the  lower 
end  if  we  use  p-1  instead  of  p  frequencies,  etc. 

For  a  given  K  and  maximum  PRF  f^,  the  highest  blind-speed  region  which 

we  wish  to  unblind  is  centered  at  (K-l)f^.  We  are  given  p  and  K  from  our 

automatic  design  procedure.  The  maximum  PRF  is  f  =  c/ (2K  )  and  we  a  I  wavs 

r  a 
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choose  this  as  our  first  PRF.  We  pick  our  second  PRF  in  the  following  marine : . 
It  is  chosen  so  that  it  blinds  the  region  just  above  the  blind  region  cen¬ 
tered  at  (K-l)f  .  Thus  we  unblind  the  region  centered  at  (K-l)f^.  This 
second  PRF  is  c^osen  to  as  large  as  possible  but  no  larger  than  f  ^ . 

The  blind  region  centered  at  (K  -2)f  is  checked  to  see  if  it  is  now  un- 

p  r 

blinded  when  f^  is  present.  If  it  is  unblinded,  we  similarly  check  (  F-  3  )  f  ^  , 
etc.,  until  we  reach  a  region  where  there  is  some  blinding,  say  at  jf  .  We 
compute  where  the  upper  end  of  the  blinded  region  lies.  The  frequency  f^  is 
chosen  so  that  it  blinds  the  region  just  above  this  upper  end.  Thus  we  un¬ 
blind  the  region  in  the  vicinity  of  jf r .  Again  f^  is  chosen  to  be  as  large 
as  possible  but  no  larger  than  f^. 

This  procedure  will  generate  decreasing  frequencies  f.  for  i  =  l  ,  ...,  p 
(where  f  =  f^).  When  a  blind  region  exists  in  the  vicinity  of  ji  the  above 
procedure  implies  we  choose  the  frequency  f^  as 


-i  fr  j  +  1 

Here  A  is  the  width  of  the  blind  zone  and  is  given  by 


(f  -  A  +  v) 


(4-1) 


2Rax(K,p) 

We  thus  use  R  and  the  computer  routine  which  computes  x(K,p)  to  obtain  A 
as  a  function  of  the  real-time  values  for  (K,p). 

The  above  procedure  assumes  that  if  there  was  blinding  in  the  vicinity  of 
jf  ,  the  choice  of  f^  in  (4-1)  totally  unbliiuled  t  lie  legion  near  jf^. 

This  total  unblinding  is  true  except  for  cases  where  both  A  is  large  and  j  is 
small.  When  these  cases  exist  we  repeat  our  unhlinding  procedure  in  the 
vicinity  of  jf  until  enough  f.  are  used  to  totally  unblind  the.  region  near 

Jfr- 

In  Figure  4-2  we  show  the  blinded  region  in  the  vicinity  of  ji  .  Thus  in 
the  expression  for  f  ^  we  have  that  v  is  the  width  of  the  unhli tided  part  at  the 
top  of  the  original  blinded  region.  Also  u  is  similarly  the  width  of  the  mi 
blinded  part  at  the  bottom  of  the  original  blinded  region.  The  region  we 
show  blinded  may  not  be  totally  blinded.  It  may  have  gaps  whirl)  uie  not 
blinded.  However,  u  and  v  define  the  region  where  there  is  some  blinding. 


A 
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We  have  generated  a  FORTRAN  program  which  determines  u  and  v  at  Jf  give', 
that  frequencies  f  ,  .  .  .  f  ^  are  present.  We  have  used  this  inutlne  and  ! 

to  determine  the  frequencies  needed  to  unblind  the  regions  oi  concern.  The 
program  can  also  be  used  to  determine  the  u  and  v  which  exist  when  we  do  not 
use  enough  frequencies  to  unblind  the  total  region  below  Kf^. 

The  FORTRAN  program  is  shown  in  the  Appendix.  The  main  subrout in<  is 
called  BLIND  (K ,  P ,  FR,DI )  .  For  inputs  of  K,p  and  FR  it  determines  t  lie  frequences 
to  be  used,  any  regions  which  remain  blinded  and  the  number  of  frequencies 
actually  used  if  less  than  p.  If  DI=0  then  the  pi  op.  ran:  computes  A  from  the 
optimum  relationship,  i.e., 


A 


FR _ 

x(K.,p) 


(4-1’) 


If  DI/O  then  it  uses  this  value  as  the  value  for  A.  The  philosophy  in 
using  (4-1 )  to  determine  the  f  is  that  we  wish  to  unblind  a  region  usiiij 
as  high  a  frequency  as  possible  and  at  the  same  time  have  it  blind  regions 
that  we  know  previous  frequencies  have  unblinded.  Out  experience  with  the 
algorithm  in  the  run  using  many  cases  is  that  at  worst  it  blinds  d'.".  oi  the 
lower  end  of  the  lowest  blind-speed  region.  None  of  the  higher  blini  ‘peed 
regions  ever  remained  blinded.  This  of  course  assumes  that  we  use  the  >  ,t  i- 
muni  relationship  for  A  in  (4-2). 

In  most  cases  that  were  run,  the  blinding  did  not  exist  at  all  or  was 
well  below  the  20%  value  in  the  lowest  blind-speed  legion.  For  many  cases , 
particularly  for  higher  p  fp  >  6)  we  did  not  need  to  use  all  the  f i equine  tes 
available  in  order  to  unblind  the  total  region  of  concern. 

4.2  PERFORMANCE  IMPROVEMENT  FROM  ADAPTIVE  TECHNiqilKS 

We  wish  to  illustrate  the  payoff  from  an  adaptive  radar  by  giving  some 
examples  of  nonhomogoneous  environments.  The  payoff  occurs  in  the  teduction 
of  the  dwell  time  required  per  search  frame.  We  consider  the  case  of  a  shoit- 
range  search  which  in  the  worst  case  is  designed  fur  a  range  of  R0  km  and 
a  clutter  width  AR  ^  of  20  m/sec.  This  cluttei  width  corresponds  to  ground 
clutter  in  the  presence  of  the  worst  expected  rain  clutter.  The  temaining 
parameters  which  are  fixed  are  given  us 
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v  -  1000  in/sec 
max 

w  =  3  re 
a 

k  =  1 
P 

k  =  1 
a 

6  =  k/2 

s 

T  =  <  2  sec 
s 


If  we  take  T 

s 


while 


1.56  sec  (leaving  .44  sec  for  track)  then  we  have 


6  = 


k  v  T 
a  max  s 

k  w  6 
pas 


331 


P 


AR  „/v 
c£  max 


.02 


For  these  values  of  6  and  p  the  subroutine  OPK  yields  the  maximum  value  of 
K-15  along  with  the  corresponding  value  of  p=6.  The  wavelength  of  trans¬ 
mission  is  then  given  by 


4v  R 

_  max  a 

cK 


7.1  cm 


Thus  we  operate  the  transmitter  at  a  wavelength  of  7.1  cm  and  in  the  worst 
case  need  a  dwell  time  for  search  of  1.56  sec. 

We  next  divide  the  search  sector  inLo  two  equal  parts  so  that  each  sector 
has  a  width 

AO  =  0  n  =  7i  /  4 
s 


and  we  perform  a  two-zone  search.  We  do  this  to  examine  how  a  nonhomogeneous 
environment  can  appreciably  reduce  the  dwell  time  required  per  search  frame. 

As  our  first  example,  suppose  R  Is  80  km  in  the  first  zone  while,  due  to 
the  change  in  terrain,  it  is  only  50  km  in  the  second  zone.  We  still  assume 
the  worst  clutter  width  of  20  m/sec  for  both  zones.  The  wavelength  of  the 
transmitter  must  stay  fixed  at  7.1  cm.  Thus  we  obtain  K  from 

4v  R 

m  ri  v  u 
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In  the  first  zone  K  is  still  15.  In  the  hckuhi  .  ...  f  'll. in.  w<  nun.: 

use  K=10  so  that  we  extend  R  to  53  km  in  th<-  set enu  tout  . 

In  the  first  zone  we  then  have  K-15,  p-.O.'  while  in  i  h.  second  />.m 
K  =  10,  p  =  .02.  We  now  use  subroutine  OPB  to  find  the  minimum  r  and  e  t  r  respond  i  np 
p  in  each  zone  for  the  given  K  and  p.  We  then  use 


T  . 
si 


(-  k  w  AG 
P  a  . 

k  v 
a  max 


1  =  1.2 


to  find  the  dwell  time  needed  for  search  in  each  azimuth  zone.  The  total 
dwell  time  is  then  given  by  T  =  T  +  T  lor  the  short-range  bursts. 

The  result  of  the  above  procedure  for  this  example  as  well  as  other  ex¬ 
amples  are  given  in  Table  4-1.  We  see  from  this  table  that  the  dwell  time 

can  be  decreased  by  anywhere  from  about  one-half  second  to  one  second  if  tin- 

nonhomogeneous  nature  of  the  environment  is  taken  into  account.  We  further 

see  that  for  the  four  cases  considered  the  maximum  of  R  and  AH  n  for  each 

a  c  x 

case  are  80  km  and  20  m/sec.  These  are  the  values  in  the  worst  case  design. 
Thus  we  have  achieved  significant  improvement  over  the  worst  c..se  even  though 
one  or  the  other  zone  has  worst  case  conditions  associated  with  it. 

We  have  given  an  example  of  dividing  a  region  into  two  zones.  Ideally 
we  should  divide  the  region  into  many  zones  so  that  in  the  limit 


where 


region 


dT 

s 

dO 


s(n;k 


w 

a 


k  v 
a  max 


Here  °(Q)  is  obtained  by  substituting  (K,p)  into  subroutine  OPB.  We  have 
that 

4v  K  (0) 

K  =  ,naX  a  -  . 
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while 

p  ■  a^«(0) 

v 

max 

As  a  practical  matter,  we  use  sums  and  zones  (rather  than 
differentials)  because  we  require  integer  mathematics  (e.g., 
of  pulses).  However,  as  seen  from  the  examples  in  Table  4-1, 
the  region  into  a  small  number  of  zones  can  yield  appreciable 
the  required  dwell  time  for  scanning. 


integrals  and 
nteger  number 
even  dividing 
improvement  it 
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5.  CONCLUSIONS 

The  major  conclusion  that  is  evident  from  the  preceding  results  is  that 
the  potential  payoff  from  adaptivity  is  great.  The  radar  can  be  designed  to 
meet  average  conditions  rather  than  worst-case  conditions,  with  the  knowledge 
that  the  radar  will  adapt  to  changing  conditions.  We  have  found  in  our  ex¬ 
amples  that  large  reductions  in  search  dwell  time  and/or  operating  wavelength 
can  result  from  the  adaptivity. 

The  major  remaining  question  is  how  much  improvement  can  be  realized  in 
practice.  Although  we  have  already  defined  the  basic  design  for  the  adaptive 
radar,  we  cannot  yet  answer  the  difficult  performance  question.  The  problem 
is  that  the  degree  of  performance  improvement  depends  on  the  range  of  environ¬ 
ments  to  be  encountered.  It  therefore  becomes  important  to  simulate  the  adap¬ 
tive  radar  under  a  broad  variety  of  conditions  in  order  to  define  performance. 

The  simulation  of  the  radar  system  is  precisely  the  next  task  that  is 
planned.  As  Illustrated  in  Table  5-1,  the  next  task  (Task  3)  is  about  eight 
months  in  duration  and  will  involve  a  testing  of  the  simulated  system  against 
realistic  environments.  We  expect  to  gain  from  this  simulation  not  only  an 
estimate  of  performance  but  also  a  feeling  for  the  measurement  accuracies  re¬ 
quired  of  the  environmental  sensor  that  measures  the  adverse  clutter  environ¬ 
ment.  The  final  task  (Task  4)  will  then  involve  the  actual  design  of  the  entire 
system  with  such  an  environmental  sensor  included. 
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Following  are  the  codes  for  the  various  subroutines  and  functions  that 
described  in  the  text. 


A- 2 

SUPKCUIIM  CPKlHtU.Hl'C.K.P.Bt  TAN  I 
IMICH  P 

CCrrLh  /TABLES/  *160, 2?) 

If  l  BhO-Cf  .1.1  THE  N 
H-C 
P-C 

HE  1  a  r- PE  1  A 
BE  llHf 
EfC  If 

BUM-  1-/ 12. M  1.-RH0I  » 

IF  iBFTA.LL.BTESI)  THEN 
K-C 
P-C 

BE  I  A  r- BE  T A 
R  f  TURN 
£  NC  IF 

IF  (RHC.CT..25.CR.BE TA.LT.6. >  THEN 
*  •  1 
P  •  1 

eil *r-  1  •/  I  2. *  c l.-RhC  > I 

B  f  TURN 
E  NC  IF 

IF  (BHC.EC.O.)  THEN 
KS-6C 
CC1C  ICC 
ENC  IF 

KS-  INT  (  1  .  /  (  2  •  •BHO)  I 
KS  -MN04K5  ,60 
20C  I-  1  •  / I RHC • *  I KS .121  I 

If  I ! •  LE  -FLOAT  IKS  I  )  Thf  N 
KS  -KS-  1 

If  IKS  .(0.  1  |  THE  N 
K-  1 
P-1 

BE  TAN  -  BE  T A 
WE  TLBN 
EL  SE 

CCTC  200 
f  NC  IF 
ENC  If 

If  (KS.EC.2J  THE  N 

CALL  CP6  12  ,  RHC  .BE  T  AN , P ) 

If  IBE  TAP*. LE. BETA)  TFEN 
K*i 

CALL  SCPK  IBE  TA  ,K,RHC,P  ,BET  AN  ) 

BE  TURN 
ELSE 
K  -  1 
P-1 

PE  1  AN  -  1 .✓! 2. ♦ I  1. -RHC I  1 
RE  TLBN 
ENC  If 
E  NC  If 
IOC  KL-2 
KH  -  N  S 

CAIL  CPH ( K S  ,RHC  ,BE TAR  ,p > 

IF  I  8E  T  A  K  .1  f  .  Pf  I  A  I  THEN 
K  •  K  S 

C  A  l  l  STPmIBE  I  A  ,  K  «  R  h()  ,P  ,Hfl  ^  I 
wt  Torn 
INC  If 

JU(J  KN  -  I  K  L  'Kf  1  /2 

r  » t  »  r  » •  i  »*  ,wi  r.».  ta-.pi 

If  |KN  .f c  .KL  .  AN L  .BE T A r .1 E  . BE  T A »  THEN 
K  •  KN 

CALL  SCPK  I  BE  T  A  ,K  ,RHO  ,P  ,  E  E  T  AN  > 

BE  TURN 
f  NC  IF 

If  I  KS  .  f  C  .  Kl  .  A  NC -BE  T  A  r  .0  T  .  RE  T  A  I  THEN 
K  -  1 

P-  1 

FF IAN- 1 ,/I2. • I  1. -RHC ) ) 
be  Turn 
E  NC  If 

IF  ( HE  I  AN  .GT . PE  I  A  I  THEN 
KH  -k  N 
CL  1C  3  C  0 
I  L  SE 

K  l  -K  N 

C  C IQ-  ICC 

I NL  If 
I  NC 
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SUBROUTINE  OPB ( K» RHO.  BETAM. P ) 

C  CALLED  BY  OPK. 

REAL  KREAL 
INTEGER  P.PP 

COMMON  /TABLES/  X(60. 12) 

KREAL=FLOAT  (  K  > 

PREAL=FLOAT ( P  > 

IF  (RHO.  EQ.  O.  )  GOTO  200 
TESTK=1.  /  (  2.  *RHO> 

IF  (RHO.  GE.  1.  .  OR.  KREAL.  GE.  TESTK)  THEN 
BETAM=0. 

P=0 

RETURN 
END  IF 

200  IF  (K.  EQ.  1 )  THEN 

BETAM=1 .  /  (  2.  *(1.  -RHO)) 

P=1 

RETURN 
END  IF 

DO  400  J=2, 12 
COMP=KREAL*RHO*X (K.  J) 

IF  (COMP.  GE.  1.  )  GOTO  400 

BETA=FLOAT(J)*KREAL*X(K, J>/(2.  Ml.  -COMP) ) 

IF  ( J.  EQ.  12)  THEN 
BETAM=BETA 
P=J 
RETURN 
END  IF 
PP=J 

BETAP=BETA 
GOTO  300 
400  CONTINUE 

BETAM= 1 .  E5 
P  =  13 

300  NN=J+1 

DO  10  I=NN,  12 

BETA=FLOAT( I )*KREAL*X(K,  I)/(2.  *(1.  -KREAL * RHO* X ( K. 

IF  (BETA.  GE.  BETAP)  GOTO  100 

PP=I 

BETAP=BETA 
10  CONTINUE 

P=PP+1 
BETAM=BETA 
RETURN 

100  BETAM=BETAP 
P=PP 
RETURN 
END 
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SUBROUTINE  SOPK ( BETA. K. R HO. P. BETAN) 

C  CALLED  BY  OPK. 

INTEGER  P.PP 
COMMON/TABLES/ X ( 60/ 12) 

IF  (P.  EG.  2.  OR.  P.  EG.  1  >  RETURN 
PP=P 

BETAP=BETAM 
N -P  -2 

DO  10  1=1. N 
J=P- 1 

TEST=FLOAT (K)*RHO*X(K.  J) 

IF  (TEST.  GE.  1.  )  THEN 
P=PP 

BETAM=BETAP 
RETURN 
END  IF 

BETAN=FLOAT( J*K)*X(K.  J)/ (2.  -2.  *TEST) 
IF  (BETAN.  LE.  BETA)  THEN 
PP=J 

BETAP=BETAN 
GOTO  10 
ELSE 
P=PP 

BETAM=BETAP 
RETURN 
END  IF 
10  CONTINUE 

P=PP 

BETAM=BETAP 

RETURN 

END 
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100 


INTEGER  FUNCTION  FK ( BETA.  RHO, PNDT, KNDT,  P ) 

INTEGER  P,  PNOT 

COMMON  /TABLES/X ( 60; 12) 

IF  (P.EQ.  PNOT)  THEN 
FK=KNOT 
RETURN 
END  IF 
N--KNOT  — 1 
DO  100  1  =  1,  N 

T  =2.  *BETA/(X ( < KNOT  +  1  —  1  ),  P>*<2.  »RH0*BETA+FL0A7 (P) )  ) 
IF  (T.  GE.  FLOAT  <KNOT-t-l-I  )  )  THEN 
FK=KN0T+1-I 
RETURN 
END  IF 
CONTINUE 
FK=1 
RETURN 
END 
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FUNCTION  FB<K, RHO, P) 

INTEGER  P 

COMMON  /TABLES/X < 60.  12) 

CRUD=FLOAT <  K ) *RHO*X ( K,  P ) 

IF  (CRUD.  GE.  1.  )  THEN 
FB=0. 

ELSE 

FB=FLOAT(K*P)*X(K, P>/<2.  *(1.  -CRUD 

END  IF 

RETURN 

END 


o  o  o  o  o 
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SUBROUTINE  TABLE 

THIS  SUBROUTINE  SETS  UP  A  TABLE  OF  VALUES  OF  THE 
FUNCTION  X X  <  K,  P )  FOR  K  =  2,  .  .  .  ,  60  AND  P  =  2,  ...  12 
THIS  TABLE  IS  THEN  PUT  INTO  THE  COMMON  BLOCK 
'TABLES'.  THIS  SUBROUTINE  MUST  BE  CAI  I  FD  BEFORE 
USING  ANY  OTHERS. 

INTEGER  P 

COMMON  /TABLES/  X<60. 12) 

DO  10  1  =  1.60 
X<I,  1)=0. 

10  CONTINUE 

DO  20  1  =  1,  12 
X  (  1.  I  )  =0. 

20  CONTINUE 

DO  40  P=2. 12 
DO  30  K=2. 60 
X  (  K,  P  )  =  X  X  ( K,  P) 

30  CONTINUE 

40  CONTINUE 

RETURN 
END 


n  r>  o 
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FUNCTION  XXtK.P) 

THIS  FUNCTION  IS  DEFINED  FOR  K=2,  .  .  .  , 60,  AND 
P=2»  .  .  .  ,  12.  THE  VALUES  OF  THIS  FUNCTIPON  ARE 
USED  IN  MOST  ALL  OF  THE  SUB-PROGRAMS 
INTEGER  P 
REAL  KK 
KK=FLOAT ( K ) 

PP=FLOAT  <  P ) 

X=KK**<  1 .  /<PP-1.  ))  +  l.  /  ( KK#*  (  1.  /<PP-1.  )  )  ) 

100  TOP  =  X**P-2.  *X**<P-1 )-KK*X+KK+l 

BOTTOM=PP#X*-*-(P-l  )-2.  *<PP- 1.  )*X**(P-2)-KK 

DEL TAX --TOP /BOTTOM 

X=X+DELTAX 

IF  (  ABS  (  DELTAX  ).  GT.  1.  E-12)  GOTO  100 

XX  =  X 

RETURN 

END 
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SUBROUT  I  NE  BLINIKK.F'fFRfDI  »M»UBL»DBL»F  ) 

IMPLICIT  DOUBLE  PRECISION ( A-H r O-Z ) 

DIMENSION  F( 1 ) ,S1B( 13) >S2A( 13) »S3B( 13) ,S3A( 13) ,UBL  ( 1  )  ? DDL (  1  ) 
INTEGER  F' 

COMMON/T ADLES/X < 60 .  12) 

CALL  XMITD(-60»0.  f  DDL  ) 

CALL  XMITD(-60f 0.  »UBL  ) 

LL  -  2 

IF  (  K  .  N  F  .  0  )  GOTO  510 
F  YF'E  GOO 

UOO  FORMAT  (IX.  7HF  1  -  O.) 

RETURN 

510  IE  (K.NE.l)  GOTO  530 

T YF'E  520  .  FRfP 

520  FORMAT  (1X.4HFR  = , D1 6 . 8 . 2X . 4HP  =  »  12 ) 

RETURN 

530  IF  (DI.EO.O.)  DI-FR/X(K.P) 

J=K  - 1 
U-0  . 

0-0. 

F( 1 )=FR 

F(2>=F0(U»0.J.F"(  1  )  r  D I  ) 

M  =  2 
L0  =  2 

IF  <F  (  1  )  .L  T  .  (DI/JT2  .  *DI  )  )  I  0-1  01 
IF  (LO.EQ.l)  GOTO  550 
IF  (J.NE.l)  GOTO  550 
CALL  CLEAR  J  (  DI  »F  .M.J.ZIJ.ZD) 

RETURN 

550  IF  (M.EG.F)  LL=1 

DO  580  L=LO  r ( K- 1 ) 

J^K-L 

560  CALL  STATES  (DI»F»M.J.S1E''.II»  S2A  .  J,J  »  S3B  .  S3A  .  KK  ) 

NN -  NS ( II . JJ.NK) 

IF  (NN.E0.8)  GOTO  580 

CALL  EX  T (SIB. 1 1 . S2A  »  JJ » S3D  *  S3A  . KK  rBiAf DCC  »  ACC ) 

CALL  UO ( BI . B . A . DCC . ACC  .  S  3  B  » 53A  . KK . NN  .  U  .  0  ) 

IF  (U.EQ.DI .OR. (UFO) .EO.DI )  GOTO  580 
IF  (LL.NE.l)  GOTO  570 
UDL ( J ) -U 
DDL( J) = D I -U-0 
GOTO  580 
570  M  =  M  f 1 

F  (  M  )  --  F  0  ( I J  t  0  » ,J  .  F  (  1  )  ,ni  ) 

IF'  (M.F  O.F  )  LL-1 

IF  (0.1  I . (DI/J42.+DI  I  (1)))  GOTO  560 
580  CONTINUE 

RE TURN 
END 


A-  1  0 


OUOKOU  I  I  N I  U  1  AR  (III  » h  rt'ir  J  G  t  Jl  r  UOL  r  HM.  / 

iwunr  oouole  precision  (A-HrO-z) 

0 I  MENS  1  UN  F  (  1  > r  UBL (l)r  UHL ( 1 ) 

TYPL  Li 00  r  HI  r  Nr  (F  (  I  )  r  ID  r  M) 

U00  FORMAT  (///,'  01  ='r016.Gr'  M  ='r!3r/r'  F  r 13(016. 8r IX 

TYPE  S 1 0 

MO  FORMAT  (//»'  J  '  r  1  OX  r  '  UOL  '  r  10X  r  '  DHL  '  r  /  r  1 X  r  29  ( 1 H- )  ) 

00  030  J--JS.JF 

CAL  L  r.Lf:  AkJ(  01  rF  r  Mr  JrZUrZO) 

UOL ( J ) - Z  U 
OHL ( J )  -ZO 

IF  (ZO.  1-0.0.)  GOTO  030 
TYPE  030  r  J  r  UOL ( J ) rOOL  ( J) 

330  FORMAT  (  1  X  r  1 3  r  3  (  5X  r  01 6 . 8  )  ) 

030  CONTINUE 
RETURN 
E.  NO 


SUKRUUT  INF.  CLEAR  J  (  01  .  F  r  M  ,  J  r  Zl)  r  ZO  ) 

IMPLICIT  00 UOL F".  PRECISION  (A-HrO-Z) 

0 l ME  NSI ON  F ( 1 ) r  S 1 0 ( 1 3 ) r  S3 A ( 1 3 ) r  S30 ( 1 3 ) r S3 A (13) 
CALL  STATES (01 rFrMr Jr  S3  hr 1  1  rSOAr JJrS30rS3ArKK ) 
NN  N5(  I  I  r JJrKK) 

IF  (NN.cn, 8)  GOTO  SOO 

CALL  EXT  ( SI  Or  T I  r S3A r J J r S30 r S3 A r KK r 0 r A r DCC  r ACC ) 
CALL  UU ( 01  r  0 , A  r OCX  >  ACC  r  S30  r S3 A r  KK  r  NN  r  U  r  V } 

IF  (U.LO.OI .OR.U+V.EQ.OI )  GOTO  SOO 
ZU  U 

ZO  O  l  (J  0 
RF TURN 
ZU  0. 

/ ; 1  0  • 

RL I  URN 
FNTl 


SOO 
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SUBROUT I NE  EXT  <  S 1 B » I » S2 A » J , S3B » S3A » K f B » A » BCC , ALT  ' 
IMPLICIT  DOUBLE  PRECISION  <A-H»0-Z) 

DIMENSION  SIB (13) r  S2A< 13) » S3A( 13 ) rS3D( 13) 

B=0 . 

A=0. 

BCC=0. 

ACC=0. 

IF  (I.EQ.O)  GOTO  510 
B=S1B( 1 ) 

IF  (I.EO.l)  GOTO  510 
DO  500  I 1=2 » I 

B=DMIN1 (B»S1B(II)) 

500  CONTINUE 

510  IF  (J.EQ.O)  GOTO  530 

A=S2A<1> 

IF  (J.EQ.l)  GOTO  530 
DO  520  JJ=2» J 

A=DMIN1(A»S2A< JJ>  > 

520  CONTINUE 

530  IF  (K.EQ.O)  RETURN 

BCC=S3B< 1 > 

ACC=S3A<1) 

IF  (K.EQ.l)  RETURN 
DO  540  KK=2»K 

BCC=DMIN1(BCC»S3B<KK) ) 

ACC=DMIN1 < ACC»S3A(KK) ) 

540  CONTINUE 

RETURN 
END 


FUNCTION  FB<NfRHO»P) 

IMPLICIT  DOUBLE  PRECISION  <A-H»0-Z) 
INTEGER  P 

COMMON  /TABLES/X < 60f 12) 

CRUD=FLOAT <  K ) *RHO*X  < K»P) 

IF  (CRUD.LT.l.)  GOTO  500 
FB=0. 

GOTO  510 

500  FB=FLOAT  <  N*P  >*X<KrP)/<2.*<l. -CRUD ) ) 

510  RETURN 

END 
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INTEGER  FUNCTION  FN < BETA » RHO »PNQT  rKNOT rP> 

IMPLICIT  DOUBLE  PRECISION  <A-H»0-Z> 

INTEGER  PtPNOT 
COMMON  /T ABLES/X  <  60 » 1 2 ) 

IF  (P.NE.PNOT)  GOTO  500 

FK=KNOT 

RETURN 

500  N=KN0T-1 

DO  510  1=1 tN 

T=2.*BETA/(X( (KNOT+l-I ) »F>#(2.*RH0*BETA+FL0AT(P> ) ) 
IF  (T.LT .FLOAT (KNOT+l-I ) )  GOTO  510 
FK=KN0T+1-I 
RETURN 
510  CONTINUE 

FK=1 
RETURN 
END 


FUNCTION  NS<ItJ,K) 

IMPLICIT  DOUBLE  PRECISION  (A-HrO-Z) 

IF  (I+J+N.NE.O)  GOTO  500 
NS=8 
RETURN 
500  N5=l 

IF  < I . GT  *  0 . AND . J . GT  « 0  *  AND ♦ K . GT  .  0 )  RETURN 

NS=2 

IF  (I.GT.O.AND.K.GT.O.AND. J.EG.O)  RETURN 
NS=3 

IF  ( K • GT  4  0 . AND ♦ J • GT • 0 • AND . I . EQ«  0 )  RETURN 
NS=4 

IF  (I.GT.O.AND. J.GT.O.AND.K.EQ.O)  RETURN 
NS=5 

IF  ( K . GT . 0 . AND . J . EQ . 0 . AND • I . EG • 0 >  RETURN 
NS=6 

IF  (I.GT.O.AND. J.EO.O. AND. K.EQ.O)  RETURN 
NS=7 

RETURN  j 

END  1 
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SUBROUT  I NE  OF'B  <  K  t  RHO p BETAM p  F*  > 

IMPLICIT  DOUBLE  PRECISION  <A-HpO-Z> 

C  CALLED  BY  OF’K. 

REAL  KREAL 

INTEGER  FpPP 

COMMON  /TABLES/  X<60p12> 

KRE AL=FLOAT  <  K ) 

F'F<EAL=FLOAT  <  P  > 

IF  (RHO.ECJ.O.)  GOTO  200 
TESTK=1 ./<2.#RH0> 

IF  (RHO.LT.l. .AND.KREAL.LT.TESTK)  GOTO  200 
BETAM=0 . 

P=0 

RETURN 

200  IF  (K.NE.l)  GOTO  201 
BET AM=1 . / ( 2 . *  < 1 . -RHO ) ) 

P-1 

RETURN 

201  DO  400  J=2 p  12 

COMP=KRE AL*RHO*X  <K» J> 

IF  <  COMP . GE . 1 . )  GOTO  400 

BETA=FLOAT  <  J  >  *KREAL*X  <  K»J>/(2.*<1. -COMP ) ) 

IF  (J.NE.12)  GOTO  202 

BETAM=BETA 

P=J 


RETURN 

202 

FF=J 

BETAP=BETA 
GOTO  300 

400 

CONTINUE 
BETAM=1 • E5 
P=13 

300 

NN= J+l 

DO  10  I=NNp 12 

BETA=FLOAT < I ) *KREAL#X (KpI>/(2.*C1. -KREAL#RHO#X ( K p I ) > ) 

IF  <  BETA . GE .  BETAF' )  GOTO  100 

FF=I 

BETAF=BETA 
10  CONTINUE 

F’-PF'+l 
BETAM=BETA 
RETURN 

100  BETAM=BETAP 

F=PP 
RETURN 
END 
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SURKOUI  INI  1)1  RtRtl  A*KIIU.K.I  .RCTAM) 

lm-Licn  ruugli:  ikecision  ca-h.o-z> 

IMTLGER  I 

COMMON  / 1 ARLES/  XlfcO*12> 

ir  <kho.lt, i.>  goto  *00 

KsrO 

(•<■0 

ret am- re? a 

RE  I  URN 

*00  ttlEST-l./<:?.*<l.-RHOJ> 

IF  < RE  1 A . 01  .1*1  ESI  I  GOTO  *05 

R=0 

I'-O 

Rr  I  AM  RCTA 
W  )URN 

*05  ir  <RHO.LI  . .25.ANR.RCTA.GE.6. >  GOTO  *10 

K- 1 
F=1 

RET AN- 1 ,/<2.»tl .-RMO)> 

RETURN 

*10  IF  (RMO.NE.O.)  GOTO  *15 

R5=60 
GOTO  100 

*15  KS=INT<1./<2.*KH0>> 

KS=M1NOIXS,60) 

*00  T=1./<RH0*X<KS«12)> 

If  (T.OT.FLOATIKS))  GOTO  *20 
KS=KS-| 

If  <RS.NE.11  GOTO  200 

K-  1 

t~i 

RtlAM  RETA 
RETURN 

*20  If  IRS.NE.2I  GOTO  100 

CALL  OFR  <  2 » RHO  •  RE  T  AM  >  f  ’  > 

IF  (RETAM.GT.RETA)  GOTO  *50 

K.=2 

CALL  SOf  "R  <  FC  TA  »  R  »  RH0>  F  r  RE  T AM  I 
RETURN 
*50  K=1 

I  =1 

Rf  TAN-RE  TA 
RE  TURN 
100  KLC2 

KM=KS 

CALL  OIR<hSfRHO»RETANrF> 

IF  <RETAM.GT.RCTA)  GOTO  300 
K-KS 

CALL  SOFKIRETA.R.RHOtFf RET AM) 

RETURN 

300  KNMKl  FRH)/2 

CAlt  UF'R  <  KN , RHO  >  RE  T  AH » P  > 

IF  <KN. NE .  KL .OR . RET AM. OT • RET A)  GOTO  *60 
K  RN 

C  A\  L  r.QF'K  <  RE  T  A  *  K  *  RHO  .  f  • »  HE  T  AH  > 

RETURN 

*60  1)  (RN.NE  .KL  .OR «  RET  AM  .IE  •  RET  A)  GOTO  *70 

R»T 
F  =  1 

RE  T  AM -t«/(2.t(l. -RHO ) ) 

Rl  TURN 

*70  IF  <RI  1AM. If  .RETA)  GU10  *80 

KM-RN 
GOTO  300 
ABO  KL  --KN 

GOTO  300 
FNR 
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« 


10 

20 

40 


50 


AO 

30 


C 


20 


30 

10 


'  I  ll/ll  I  HIM  KMX.  T.h.A.M-  > 

IMPLICIT  RIJURU  II.IL1MON  (A-H.O  /> 
I'lMLN'JICJN  XUJl.VUJl 
1*0  10  1-1  tKK 

If  (YU  >  .HI  .1*)  GOTO  1*0 
CONI INUf 
KR=  H 
RETURN 

IF  (XU  )  .LC.Rl  0010  40 

RR=R 

F.ETURN 

IF  (Y(l).LT.AI  GOTO  50 
KR=A 
RE TURN 
KM  YU  ) 

IF  CKK.EO.  1  .OR.  1  .EC1.KK)  RETURN 

I  TEMP*  HI 

HO  30  K=1T£MP.KK 

IF  (X(K).GT.KI<>  KFTUKN 

IF  (YCK).LE.RR)  GOTO  30 

IF  (YCK1.LT .A)  GOTO  AO 

RR»A 

RETURN 

RR*Y<K> 

CONTINUE 

RETURN 

ENli 


SURROUT I NE  SOPK  (  RET  A  .  K  .  RHO  .  F* .  RE  T  AM  > 
IMPLICIT  HOURLE  PRECISION  (A-H.O-Z) 
CALLER  RY  OPK. 

INTEGER  F.PP 
COrtMON/T  ARLES/X ( 60  > 1 2  > 

IF  (P.t0.2.0R.P.En.l)  RETURN 
FF«P 

RETAF-RETAN 

N-P-2 

RO  10  1  =  1. N 
J=P-I 

TEST=F10ATCK>*RH0»XCK. J> 

IF  (TEST. LT. 1.1  GOTO  20 
F=PP 

RETAH=RETAP 

RETURN 

RETAN=FLOATC J*K >*X(K. J1/C2. -2. *TCST ) 

IF  CRETAN. GT.RETAl  G010  30 

PP-J 

RETAP*RETAN 
GOTO  10 
F=PP 

RETAM=HETAP 

RETURN 

CONTINUE 

P*PP 

RE TAM* RE TAP 

RETURN 

END 
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SUIU-UUI  INI  SCIK1  <  b  .  »■  .UI*.UI.«I»I.«I'I  > 
IMPLICIT  IHJUKLf  1-M.riblUM  <A  H.0-2) 
I‘J  NLNSJON 

lio  10  J=1  #KK 
K<  J> 

\0  CONI 1NUE 

OI< » 1  >  *  !>  *  I  > 

UA<  1  >»M  1  ) 

IF  IKR.EO.l)  RETURN 
JiO  bO  1=2,RK 

BO  20  K*1  »  < 1-1 1 

IF  IS(I).LI.UMM)  GOTO  30 
20  CONTINUE 

uMi^sm 

UA< I >=M  I ) 

GOTO  00 

30  I>0  *0  L=K,  <1-1  > 

UH<L+1)=UR<L> 

UA<L-*1)=UA<L) 

40  CONTINUE 

UR <  K ) =S< I > 

UA<K>=M  J  > 
bO  CONT 1 NUE 

RE  T  URN 
ENH 


10 

20 

30 


SUBROUTINE  STATES!  t*I»r  .  M.  JJ.Sl  B.  I  r  S2A .  J,  S3R .  S3A.K  ) 
IMPLICIT  BOURLC  PKECIS1DN<A-H.0-Z> 

MMENSION  F<  J  > . SI R< 1 3 > . S2A C 1 3 > . S3R <  1  3 >  ,S3A(13) 

1=0 

J=0 

K*0 

HO  100  11=2, M 

0=FL0AT< JJ)*F<1»/F< JJj 
N=  INI  <01 

BC  =  FLOAT  <  N)  *F<  1 1  MBJ 
AC  =  RC+F<Il>-2*Ii| 

TEMP=FL0AT<JJ)*F<1) 

IF  <»C.GT.1EHP.0R.AC.LT.TENF>  GOTO  10 

1=0 

J=0 

K=0 

RETURN 

IF  <BC.LE.TEMP.OF. AC. GE. TEND  GOTO  20 
K  =  R  +  1 

S3H<R')  =  RC-TEMP 
S3A<K1=-AC+TEMP 
GOTO  100 

IF  (BE. LE. TEMP)  GOTO  30 
1  =  1  +  1 

S3  B< I ) »BC-T CMP 
GOTO  100 
J=  J+l 

G2A<-IV.  -ACT1LMF 
CONTINUE 
M  I  URN 
LNH 


100 
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SUBROUTINE  TABLE 

IMPLICIT  DOUBLE  PRECISION  (A-HfO-Z) 

C  THIS  SUBROUTINE  SETS  UP  A  TABLE  OF  VALUES  OF  THE 

C  FUNCTION  XX<KfP>  FOR  K=2f...f60  AND  P=2f...f12. 

C  THIS  TABLE  IS  THEN  PUT  INTO  THE  COMMON  BLOCK 

C  ’TABLES*.  THIS  SUBROUTINE  MUST  BE  CALLED  BEFORE 

C  USING  ANY  OTHERS. 

INTEGER  P 

COMMON  /TABLES/  X(60f12> 

DO  10  1*1 t 60 
X(If1)=0. 

10  CONTINUE 

DO  20  1=1 f 12 

X ( 1 F I ) =0. 

20  CONTINUE 

DO  40  F‘=2  f  12 
DO  30  K=2f60 
X(KfF‘)=XX(KfP) 

30  CONTINUE 

40  CONTINUE 

RETURN 
END 


SUBROUTINE  UV<DI fBfAfBCCfACCfS3BfS3AfKfNfUfV) 

IMPLICIT  DOUBLE  PRECISION  <A-HfO-Z> 

DIMENSION  S3A< 13) fS3B( 13) 

IF  (N.GT.7.0R.N.LT. 1 )  GOTO  80 
GOTO  <10f20f30f40f50f60f70>  N 
10  CALL  UV1<DIfBfAfBCCfACCfS3BfS3AfKfUfV) 

RETURN 

20  CALL  UV2<DIfBfBCCfACCfS3BfS3AfKfUfV> 

RETURN 

30  CALL  UV3(DIfAfBCCfACCfS3BfS3AfKfUfV) 

RETURN 

40  CALL  UV4<DI fBfAfUfV) 

RETURN 
50  U=0. 

V=0. 

RETURN 
60  U=0. 

V=DI-B 
RETURN 
70  U=DI-A 

v=o . 

RETURN 

80  TYPE  90 f  N 

90  FORMAT  <1Xf37HN  OUT  OF  BOUNDS  IN  SUBROUTINE  UV.  N  =fI3> 

STOP 
END 


SUBROUTINE  UV1  <  DI »  B  >  A » BCC.  »  ACC  » S3B  r  S3A » K r U »  V  > 
IMPLICIT  DOUBLE  PRECISION  <A-HfO~Z> 
DIMENSION  S3B< 13) t S3AC13) ,UB< 13) »UA< 13) 
DAC=DI-ACC 
D»C=DI-BCC 

IF  < <A+B> .GT.DI)  GOTO  10 

U=PI 

V=0. 

RETURN 

IF  < B . LE ♦ DAC . OR . A . LE . DBC >  GOTO  40 

U=DI-A 

V--DI-B 

RETURN 

IF  <  B . GT . DAC . OR . A . GT . DBC )  GOTO  40 
CALL  SORT (S3BrS3A»UB>UArKrDI) 

U^RB ( UB  t UA  >  D I  -  A » B  »  K  > 

IF  (U.NE.B)  GOTO  30 

U=DI 

V=0. 

RETURN 

CALL  SORT <S3Af S3B»UBf UA»Kr  EH  > 
V=RB(UB,UA>DI-B>DI-U»K) 

RETURN 

IF  <B.LE.DAC.OR.A.GT.PBC>  GOTO  SO 
CALL  SORT <S3B*S3A*UBrUA*KrDI > 

U=RB ( UB , UA  r  DI -A  » DAC  »  K ) 

V=DI-B 

RETURN 

CALL  SORT (S3ArS3BrUBrUArKrDI > 

V=RB  ( UB  f  UA  r  DI  -B  »  DBC  »  K ) 

U=DI-A 

RETURN 

END 


SUBROUTINE  UV2 < HI r B > BCC t ACC  * S3B * S3A r K t U *  V ) 
IMPLICIT  DOUBLE  PRECISION  (A-HrO-Z) 

DIMENSION  S3B< 13) » S3A  ( 13) r UB ( 13 ) » UA < 13) 

U=0 . 

DAC-DI-ACC 

IF  <  < B ♦ GT . DAC . OR . B . GT . BCC ) • AND « B . LE . DAC )  GOTO  10 

V=DI-B 

RETURN 

CALL  SORT  <  S3A » S3B  »  UB » UA  r  N  *  I'l  > 

V=RB ( UB  t UA  f  D I -B  t DI - BCC  t K ) 

RETURN 


non 
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SUBROUTINE  UV3 <  HI f  A »  BCC » ACC  t S3B f S3A , K , U »  V ) 
IMPLICIT  DOUBLE  PRECISION  <A~HfU  Z> 

DIMENSION  S3B <  13  >  f S3A <  1 3 )  f UB ( 13 ) »UA( 13) 

V=0. 

DBC=DI-BCC 

IF  < (A.GT.DBC.OR.A.GT.ACC) . AND . A . LE . DBC >  GOTO  10 

U=DI -A 

RETURN 

10  CALL  SORT  < S3B f S3Af UB fUA f K f DI ) 

U=RB (UB f UA f  DI-A f DI-ACC  f K ) 

RETURN 

END 


SUBROUTINE  UV4 (DIfBfAfUfV) 

IMPLICIT  DOUBLE  PRECISION  <A-HfO-Z> 

IF  ( (B+A) .GT.DI)  GOTO  20 

U=DI 

v=o. 

RETURN 
20  U=DI-A 

V=DI-B 
RETURN 
END 


FUNCTION  XX(KfP) 

IMPLICIT  DOUBLE  PRECISION  (A-HfO-Z) 

REAL  KKfPP 

THIS  FUNCTION  IS  DEFINED  FOR  K=2f...f60f  AND 
P=2 f . ♦ • f 12.  THE  VALUES  OF  THIS  FUNCTIPON  ARE 
USED  IN  MOST  ALL  OF  THE  SUB-PROGRAMS. 

INTEGER  P 
KK=FLOAT  <K) 

PF'=FLOAT  (  P ) 

X=NN**  < 1 . / < PP- 1 . ) )  +  1 . / <  NK**  < 1 . / <  PP-1 . ) ) > 

100  T0P=X**P-2.*X**<P-1)-NN*X+KK+1 . 

BOTTOM=PP*X**  <  P-1 ) -2 . *  <  PP- 1 . ) *X** ( P-2 ) -KK 

DELTAX=-TOP/BOTTOM 

X=X+DELTAX 

IF  (ABS(DELTAX) .GT.l.E-12)  GOTO  100 

XX=X 

RETURN 

END 
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